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EXECUTIVE  SUMMARY 


Studies  of  processes  that  determine  the  ion  compositions  of  the  iono¬ 
sphere  are  of  great  practical  interest  in  the  field  of  long-range  optical  and 
radio  communications.  Apart  from  the  major  positive  ions  N0+  and  O2, 
relatively  large  amounts  of  metal  ions  are  found  in  the  D-region.  It  is 
believed  that  metal  atoms  are  deposited  by  ablation  from  meteors  and  ionized 
by  photoionization  and  charge  transfer  reactions. 

Emission  from  neutral  and  ionic  forms  of  metals  have  been  identified  in 
meteor  spectra  for  sodium  (Na),  potassium  (K),  lithium  (Li),  magnesium  (Mg), 
iron  (Fe),  and  calcium  (Ca).  In  addition,  ions  of  Fe,  Mg,  Na,  and  Ca  have 
been  detected  using  a  rocket-borne  ion  mass  spectrometer. 

Quantitative  results  on  Na  emissions  from  ground-based  measurements  have 
documented  the  seasonal  and  diurnal  variations  of  the  neutral  sodium  layer 
around  90  km.  Among  the  metals,  sodium  is  the  most  studied  because  of  the 
large  amount  of  data  available  on  both  neutral  and  ion  distributions.  Recent 
efforts  at  modeling  the  metal  ions  in  the  D-region  have  been  helped  by 
incorporation  of  a  relatively  sophisticated  positive  and  negative  ion 
chemistry  in  the  calculation  of  the  concentrations  of  the  major  positive  and 
negative  ions.  However,  attempts  at  studying  the  detailed  diurnal  behavior  of 
the  ions  are  hampered  by  the  lack  of  knowledge  of  the  diurnal  variations  of 
the  neutral  species  such  as  0,  O3 ,  and  NO. 

A  recent  suggestion  that  certain  atmospheric  metals  (particularly  sodium) 
may  play  an  important  role  in  the  stratospheric  chlorine  chemistry  casts  the 
studies  of  ion/neutral  chemistry  of  metal  in  a  new  light.  Below  90  km,  the 
bulk  of  the  metal  species  are  in  their  neutral  form.  Once  in  the 
stratosphere,  the  metal  may  undergo  a  series  of  homogeneous  and  heterogeneous 
reactions  to  form  metallic  compounds  and  are  ultimately  removed  from  the 
atmosphere  by  rain-out  or  deposition  in  the  troposphere.  The  question  remains 
whether  the  metals  of  high  reactivity  are  present  in  sufficient  abundance  as 
to  affect  the  other  chemical  cycles  in  the  stratosphere  and  mesosphere. 

Three  years  ago,  we  proposed  to  initiate  a  modeling  study  program  to  sej: 
up  a  one-dimensional  diurnal  model  of  the  stratosphere  and  mesosphere  to  study 
the  chemistry  of  the  metal  neutrals  and  ions  in  the  D  and  E  regions.  A 
necessary  first  step  in  the  project  is  to  use  the  model  to  simulate  the 


behavior  of  Che  background  trace  gases  and  ions.  Since  sodium  is  Che  species 
with  most  kinetic  data  and  observation,  we  proposed  to  concentrate  on  sodium 
using  the  available  observations  as  a  way  to  validate  the  model  results. 

Tiie  results  from  our  efforts  are  reported  here.  The  report  is  organized 
into  five  parts.  The  key  findings  from  each  are  summarized  below: 

I.  We  have  set  up  a  one-dimensional  model  for  the  stratosphere  and 
mesosphere  encompassing  the  altitude  range  60-120  km  that  calculates  the 
altitude  profiles  and  diurnal  behavior  of  the  background  trace  gases  as 
well  as  metal  neutrals  and  ions.  A  study  was  made  on  the  simulation  of 
the  altitude  distribution  and  diurnal  behavior  of  neutral  sodium  in  the 
atmosphere.  The  results  were  compared  with  the  available  observations 
with  the  purpose  of  identifying  the  photochemical  reactions  that  control 
the  behavior.  At  the  start  of  the  program,  many  of  the  key  reaction  rate 
constants  had  not  been  measured  in  the  laboratory.  Some  of  the  reactions 
that  we  discussed  have  now  been  measured.  With  these  rates,  the  calcu¬ 
lated  behavior  of  the  sodium  species  is  consistent  with  the  available  ob¬ 
servations,  indicating  that  the  assumptions  in  the  model  concerning  atmo¬ 
spheric  inputs,  and  photochemical  transformation  are  valid.  We  found 
that  an  estimated  input  of  2  x  lO1*  cm”2  sec”1  (1.2  x  102  tons/year)  of 
sodium  atom  from  ablation  of  meteor  will  produce  a  layer  of  neutral  sodi¬ 
um  atom  around  90  km  with  concentration  of  3-5  x  103  cm”3,  in  agreement 
with  observations.  Between  60-80  km,  the  sodium  species  are  in  the  form 
of  Na02  and  NaOH.  A  value  of  10_i0  cm3  sec”1  for  the  reaction  of  NaO 
with  H20  will  yield  a  diurnal  variation  of  14  percent  of  the  mean,  con¬ 
sistent  with  the  upper  bound  from  observations. 

II.  To  estimate  the  impact  of  a  metal  on  stratospheric  chemistry,  we  must 
first  establish  how  much  metal  enters  into  the  stratosphere.  The  valida¬ 
tion  of  the  sodium  results  in  the  mesosphere  support  our  estimates. 

Recent  measurements  of  the  rates  tor  reactions  of  Na02  and  NaOH  with  HC1 
and  of  the  photoabsorption  cross  section  of  NaCl  allow  for  realistic 
estimates  of  the  impact  of  sodium  chemistry  on  chlorine  partitioning  in 
the  upper  stratosphere.  Calculated  CIO  concentrations  near  50  km  could 
increase  by  up  to  a  factor  of  two  with  a  corresponding  decrease  in  calcu¬ 
lated  O3  abundances  of  as  much  as  13  percent.  The  impact  of  CIO  and  O3 

i  v 


depends  cruciaLLy  on  the  assumed  meteorttlc  source  of  sodium,  on  trans¬ 
port  from  the  mesosphere,  and  on  transformation  and  removal  of  gaseous 
sodium  in  the  stratosphere.  Since  the  fluctuations  in  the  abundance  of 
gaseous  sodium  could  induce  O3  variations  comparable  to  those  caused  by 
changes  in  temperature  and  solar  UV,  consideration  of  the  chemistry  of 
sodium  and  other  metals  may  be  an  important  element  in  efforts  toward 
early  detection  of  trends  in  upper  stratospheric  O3.  Future  increases  in 
stratospheric  chlorine  will  magnify  the  impact  of  sodium  species  on  ozone. 

III.  The  diurnal  photochemical  model  described  in  Part  I  has  been  utilized  to 

calculate  the  diurnal  variation  of  background  and  sodium  ions  between  60 
and  120  km.  Background  ion  species  include  the  primary  constituents  N0+ 
and  of*  Densities  of  Na+  are  of  order  2-3  x  102  cm-2  at  100  km,  with  a 
diurnal  variation  of  factors  of  400  and  3  at  90  and  100  km, 

respectively.  The  abundance  of  Na+  above  90  km  determines  the  topside 

scale  height  of  the  neutral  sodium  layer.  This  abundance  is  controlled 
by  the  rate  for  the  clustering  reaction  of  Na+  with  N2>  which  is  the  rate 
limiting  step  for  conversion  of  Na+  to  Na.  Downward  transport  of  Na+  due 
to  ion  drift  is  also  effective  in  lowering  Na  densities  above  90  km.  The 
model  can  successfully  simulate  the  observed  topside  scale  height  of  the 

sodium  layer  by  either  using  a  slow  rate  for  the  clustering  reaction  ~  7 

x  lO-22  cm6  sec”1  at  200°  K  or  adopting  a  downward  ion  drift  velocity  of 
20  cm/s. 

IV.  Seasonal  variations  have  been  observed  in  the  total  column  abundance  of 
Na  and  the  D-line  emission  rate.  Previous  studies  by  Garcia  and  Solomon 
(1985)  have  successfully  explained  seasonal  variations  in  0  and  O3  in 
terms  of  changes  in  the  vertical  diffusion  coefficient  Kzz  induced  by 
corresponding  variations  in  breaking  of  gravity  waves  in  the  mesosphere. 

We  have  incorporated  transport  effects  Into  our  model  and  investigated 
the  sensitivity  of  our  results  to  changes  in  Kzz  and  background  atmo¬ 
spheric  density  and  temperature.  Seasonal  changes  in  the  background 
atmosphere  induce  a  winter  maximum  in  the  Na  column  density  a  factor  of 
1.4  higher  than  the  summer  minimum  in  qualitative  agreement  with  observa¬ 
tions.  Seasonal  changes  In  Kzz  induce  a  seasonal  variation  of  tactors  of 
three  to  four  in  the  Na  column  abundance  and  the  D-line  emission,  peaking 
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at  the  equinoxes.  This  result  is  in  agreement  with  Low-latitude  D-line 
observations,  but  in  disagreement  with  high-latitude  D-line  observations. 
We  cannot  thus  reconcile  the  Na  observations  with  the  seasonal  changes  in 


Kzz  within  the  context  of  the  present  photochemical  model.  The  calcu¬ 
lated  diurnal  varation  of  the  Na  column  density  is  also  sensitive  to 
seasonal  changes  in  background  trace  species,  particularly  0  and  H20. 

V.  The  diurnal  variation  and  altitude  profiles  for  Li  and  K.  have  been  calcu¬ 
lated  using  our  photochemical  model.  Lack  of  kinetic  data  precludes 
detaiLed  analysis  of  the  behavior  of  the  above  species.  Our  results  are 
consistent  with  the  few  available  observations,  if  we  adopt  chemical 
schemes  similar  to  those  for  Na  and  assume  sources  equal  to  1U~ 1  and  10” 3 
times  that  ot  Na. 
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PART  I 


The  Diurnal  Variation 
of  Neutral  Sodium  Species 
in  the  Upper  Atmosphere 


ABSTRACT 


We  have  set  up  a  one-dimensional  model  for  the  stratosphere  and  meso¬ 
sphere  encompassing  the  altitude  range  60-120  km  that  calculates  the  altitude 
profiles  and  diurnal  behavior  of  the  background  trace  gases  as  well  as  metal 
neutrals  and  ions.  A  study  was  made  on  the  simulation  of  the  altitude  distri¬ 
bution  and  diurnal  behavior  of  neutral  sodium  in  the  atmosphere.  The  results 
were  compared  with  the  available  observations  with  the  purpose  of  identifying 
the  photochemical  reactions  that  control  the  behavior.  At  the  start  of  the 
program,  many  of  the  key  reaction  rate  constants  had  not  been  measured  in  the 
laboratory.  By  now,  some  of  the  reactions  that  we  discussed  have  been  mea¬ 
sured.  With  these  rates,  the  calculated  behavior  of  the  sodium  species  are 
consistent  with  the  available  observations,  indicating  that  the  assumptions  in 
the  model  concerning  atmospheric  inputs,  photochemical  transf orma t ion  are 
valid.  We  found  that  an  estimated  input  of  2  x  104  cm  2  sec_1  (1.2  x  102 
tons/year)  of  sodium  atom  from  ablation  of  meteor  will  produce  a  layer  of 
neutral  sodium  atom  around  90  km  with  concentration  of  3-5  x  103  cm  3,  in 
agreement  with  observations.  Between  60-80  km,  the  sodium  species  are  in  the 
form  of  Na02  and  NaOH.  A  value  of  10  ^  cm3  sec  1  for  the  reaction  of  NaO 
with  H2O  will  yield  a  diurnal  variation  of  14  percent  of  the  mean,  consistent 
with  the  upper  bound  from  observations. 
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I .  INTRODUCTION 


The  existence  ot  trace  amounts  of  atomic  sodium  (Na)  in  the  Karth's  upper 
atmosphere  was  deduced  in  the  earlier  part  of  this  century  from  observations 
of  nighttime  airglow  in  the  D  Lines  of  Na  (Chapman,  1939,  and  references  there 
given).  Subsequent  measurements  of  resonance  radiation  in  the  0  lines  were 
carried  out  at  twilight  and  during  the  day  by  numerous  investigators  in  the 
period  prior  to  1970  (see  review  by  Kvifte,  1973).  These  measurements  deduced 
the  concentration  of  atomic  sodium  in  a  narrow  layer  centered  at  around  90  xm. 
The  advent  of  laser  radar  techniques  (lidar)  has  produced  a  wealth  of  infor¬ 
mation  on  the  altitude  distribution  and  the  diurnal  and  seasonal  variation  of 
the  sodium  layer  during  the  past  L3  years  (Gibson  and  Sandford,  1971,  1972; 
Megie  and  Blamont,  1977;  Simonich  et  al.,  1979;  Clemesha  et  al.,  1982;  Granier 
and  Megie,  1982;  Gardner  et  al.,  1986).  Although  different  mechanisms  were 
originally  proposed  to  explain  the  origin  of  this  sodium  layer  (Kvifte,  1973), 
experimental  evidence  and  present  theoretical  models  favor  a  source  from  abla¬ 
tion  of  meteorites  occuring  between  80  and  100  km  (Megie  and  Blamont,  1977; 
Gadsden,  1968;  Hunten  (it  aL.,  1980). 

We  present  in  this  report  the  results  of  our  calculation  of  the  diurnal 
variation  of  neutral  sodium  compounds.  Comparison  of  our  results  with  the 
observed  properties  that  .are  determined  by  photochemica 1  processes  can  place 
additional  constraints  on  some  of  the  reaction  rates.  In  particular,  we  will 
see  that  the  reaction  of  Na()2  +  0  must  be  slow  and  that  the  lower  bound  for 
the  measured  rate  of  NaO  +  H2 < )  is  consistent  with  present  observations.  Our 
study  will  also  isolate  the  chemical  cycles  controlling  the  distribution  and 
diurnal  behavior  of  the  sodium  layer,  thus  estimating  the  sensitivity  of  our 
results  to  unknown  kinetic  rates  and  properties  of  the  background  atmosphere. 

Sections  2  and  1  summarize  the  information  about  laboratory  kinetic 
measurements  and  observations  ot  atmospheric  sodium  available  at  present, 
respectively.  Our  model  is  described  in  Section  4,  and  the  results  of  our 
calculations  are  presented  in  Section  3. 

2.  SUMMAKY  OK  KHOTOCMKMICAL  SCHKMK 

[tie  tirst  photochemi  c.a  1  model  for  sodium  compounds  was  proposed  in  1939 
by  Chapman  to  explain  the  nighttime  I)— lint*  emission.  This  mechanist,  proceeds 
through  tiie  ox  i  da  t  i  on— reduct  i  on  cycle  involving  Na  and  sodium  monoxi  (NaO) 


Na  +  O3  +  NaO  +  0 


( RR1  )* 


NaO  +  0  +  Na(2P,2S)  +  02  ,  (RR2) 

where  the  dissociation  energy  of  02  is  partially  channelled  into  the  excita¬ 
tion  of  Na  to  the  2P  state,  with  subsequent  radiation  in  the  D  lines.  The 
validity  of  the  above  mechanism  has  been  confirmed  by  kinetic  studies  and  at¬ 
mospheric  observations  during  the  past  40  years.  The  evidence  against  the 
possibility  of  other  mechanisms  has  been  recently  summarized  by  Bates  and  0 j ha 
(1980). 

The  existence  and  chemistry  of  sodium  compounds  other  than  atomic  sodium 
is  more  uncertain,  because  of  a  lack  of  laboratory  kinetic  data  and  atmospher¬ 
ic  measurements.  The  nonproton-hydrate  positive  ions  observed  in  the  middle 
atmosphere  by  Arnold  et  al.  (1978)  were  originally  interpreted  to  be  of  the 
form  Na0H2+(Na0H)m(H20)n  (Ferguson,  1978).  This  interpretation  led  Liu  and 
Reid  (1979)  to  postulate  the  production  of  NaOH  through  reactions  such  as 

Na02  +  OH  +  NaOH  +  02 ,  (RR10) 

Na02  +  H  NaOH  +  0.  (RR1L) 

where  Na02  is  produced  by  the  reaction 

Na  +  02  +  M  *  Na02  +  M.  (RR3 ) 

This  model  predicted  rapid  conversion  of  Na  to  NaOH  below  an  altitude  of 
90  km.  Direct  production  of  NaOH  from  NaO  was  proposed  by  Murad  and  Swider 
(1979)  through  the  reaction 

NaO  +  H20  -*•  NaOH  +  OH.  (RR7) 

A  value  of  (2  ±  1)  *  10  10  cm3  sec  1  for  the  rate  of  (RR7)  has  been  recently 
obtained  by  Ager  and  Howard  (1985).  Although  more  recent  measurements  of 
stratospheric  positive  ions  have  eliminated  the  possibility  of  a  sodium  core 
(Arnold  et  al.,  1981),  the  NaOH  sources  represented  by  reactions  (RR10), 
(RRll),  and  particularly  (RR7)  are  stil'  realistic. 


*RR  numbers  refer  to  reaction  numbers  in  Table  !. 


The  abundance  of  NaOH  is  also  of  interest  for  stratospheric  ozone  chemis¬ 
try.  Murad  and  Swider  (1979)  and  Murad,  Swider,  and  Benson  (1981)  pointed  out 
that  the  reaction 

NaOH  +  HC1  -*•  NaC  1  +  H20  (1) 

could  constitute  an  important  reaction  pathway  for  stratospheric  HC1  if  its 
rate  were  of  the  order  of  10  ^2  cm3  s  and  if  the  NaOH  concentration  were 
~  10-  cm  3.  Furthermore,  if  NaCl  were  sufficiently  stable  to  act  as  a  nucle¬ 
ating  agent  and  be  incorporated  into  aerosols,  reaction  (1)  could  provide  an 
important  sink  for  stratospheric  chlorine.  The  value  of  the  rate  for  reaction 
(1)  has  recently  been  measured  to  be  2.5  x  10  10  cm3  s  1  (Silver,  Stanton, 
Zahniser,  and  Kolb,  1984a).  A  more  detailed  discussion  of  the  coupling  to 
stratospheric  chemistry  will  be  given  in  Part  II  of  this  report. 

More  recently,  the  rate  for  (RR3)  has  been  measured  (Husain  and  Plane, 
1982;  Silver,  Zahniser,  Stanton,  and  Kolb,  1984b)  to  have  a  value  three  orders 
of  magnitude  higher  than  that  used  in  Liu  and  Reid's  model.  Such  high  values 
for  the  rate  of  (RR3)  raised  the  possibility  of  the  existence  of  large  amounts 
of  Na02  in  the  upper  atmosphere,  with  consequently  low  concentrations  of 
atomic  sodium.  Analysis  of  the  latest  therraochemical  data  led  Sze,  Ko,  Murad, 
and  Swider  (1982)  to  propose  the  reaction 

Na02  +  0  *  NaO  +  02  (RR6 ) 

as  a  pote nt ia 1 1 y-f as t  recycling  mechanism  for  Na02. 

Our  present  understanding  of  the  chemistry  of  sodium  compounds  is  summar¬ 
ized  in  Figure  l,  based  on  the  scheme  of  Sze,  Ko,  Murad,  and  Swider  (1982). 
Because  of  the  difficulty  in  manipulating  gas-phase  Na  in  the  laboratory,  kin¬ 
etic  data  for  reactions  involving  sodium  have  become  available  onlv  recently. 
Measurement s  are  now  available  for  the  rates  of  (RR1),  (RR3),  (RR4),  and 
c  K  K  7 ) ,  as  well  as  other  secondary  react  i  mss  involving  ii2  and  O3.  Constraints 
can  be  placed  on  th  ■  Na(-’  P )  branch  nt  (RR2  .I  i  i  u«  t  he  measured  background  con¬ 
centrations  of  O  and  <)•}  ,  the  observed  ii-line  emission,  and  kinetic  considera¬ 
tions  1 Bates  and  Ojha,  1980)  .  Addition. ii  c<.r.s  t  rn  i  nt  s  on  other  reaction  rates 
'an  be  outlined  by  comparing  ca  1  on  1  a  t  L  oils  w.!:i  1  bse  r  v  1 1  i  ms  ot  the  iistribu- 


tion  (Liu  and  Reid,  1979;  Sze,  Ko,  Murad,  and  Swider,  1982;  Thomas,  Lsherwood, 
and  Bowman,  1983)  and  diurnal  variation  of  sodium  (Kirchoff  and  Cleraesha, 

1983;  Kirchoff,  1983).  Estimated  rates  are  givr.t  in  Table  1  where  laboratory 
measurements  do  not  exist. 

In  Figure  1,  the  Na-NaO  on  the  left-hand  side  represents  the  original 
Chapman  scheme.  The  rates  for  (RR1)  and  (RR2)  are  sufficiently  fast  that  the 
[Na]/[NaO]  should  be  proportional  to  the  [0]/ [O3 ]  at  any  time  of  the  day. 

Thus,  in  the  absence  of  other  sodium  species,  the  diurnal  behavior  of  the  Na 
and  NaO  would  be  determined  entirely  by  the  diurnal  variations  of  0  and  O3. 

The  species  NaC>2  and  NaOH  act  as  temporary  reservoirs  for  the  sodium  species 
during  the  diurnal  cycle.  With  the  short  photolysis  lifetime  for  NaOH,  diur¬ 
nal  variation  of  Na  Is  expected  if  the  density  of  NaOH  is  sufficiently  high. 
The  effect  of  Na02  on  the  diurnal  variation  of  Na  would  depend  on  the  magni¬ 
tude  of  the  recycling  reactions  (RR5)  and  (RR6). 

3.  SUMMARY  OF  OBSERVATIONS  ON  NEUTRAL  IONS 

An  extensive  data  base  now  exists  providing  detailed  information  on  the 
diurnal  and  seasonal  behavior  of  the  sodium  layer.  Some  of  these  properties 
are  directly  attributable  to  transport  phenomena  on  a  seasonal  (Gibson  and 
Sandford,  1971;  Megie  et  al.,  1978;  Kirchoff  et  al.,  1981;  Kirchoff  and  Cle¬ 
raesha,  1983)  or  diurnal  time  scale  (Megie  and  Blamont,  1977;  Clemesha  et  al., 
1978;  Kirchoff  et  al.,  1981;  Clemesha  et  al.,  1982).  Nevertheless,  we  can 
deduce  from  the  data  certain  characteristics  that  have  direct  relevance  to 
photochemical  processes.  The  experimental  data  to  be  used  in  constraining  our 
model  are  as  follows: 

a)  Total  Abundance  of  Atomic  Sodium:  Lidar  observations  have  con¬ 
strained  the  diurnal  variation  of  the  Na  column  density  to  less  than 
15%  (Gibson  and  Sandford,  1972;  Clemesha  et  al.,  1982).  The  abso¬ 
lute  magnitude  of  the  total  abundance  ranges  from  a  summer  minimum 
of  ~  2  x  109  cm  2  to  a  winter  maximum  of  6-8  x  109  cm  2,  yielding  an 
average  column  density  of  A  x  109  cm  2. 

b)  Altitude  Profiles  of  Atomic  Sodium:  Lidar  observations  of  the  sodi¬ 
um  layer  at  different  latitudes  have  consistently  deduced  a  peak, 
density  of  3-5  x  iO3  cm  3  at  90-92  km  during  periods  of  low  meteor- 


itic  activity  (Megie  and  Blamont,  1977;  Simonich  et  al.,  1979).  The 
peak  height  and  density  correlates  with  the  seasonal  change  in  total 
abundance,  with  high  column  densities  corresponding  to  low  peak  al¬ 
titudes  and  high  peak  densities  (Gibson  and  Sandford,  1971). 

Diurnal  variations  in  the  peak  height  have  also  been  observed,  but 
they  are  probably  due  to  transport  effects  of  gravity  waves  (Cle- 
mesha  et  al.,  1982). 


c)  Diurnal  Variation  of  Local  Densities  of  Atomic  Sodium:  The  densi¬ 


ties  measured  by  lidar  techniques  exhibit  a  diurnal  variation  below 
84  km  (Simonich  et  al.,  1979;  Clemesha  et  al.,  1982).  The  nocturnal 


decrease  in  sodium  densities  observed  by  these  authors  is  about  a 


factor  of  2  at  83  km  and  10  at  81  km. 


d)  D-line  Emission  Rates:  Measurements  of  the  nighttime  D-line  emis¬ 
sion  in  the  southern  hemisphere  (low  latitude)  yield  minimum  values 
of  30  Rayleigh  (R)  during  winter  and  60  to  100  R  at  equinox  (Kir- 
choff  et  al.,  1981).  Mean  intensities  in  the  northern  hemisphere 
range  from  35  to  100  R  (Fukuyama,  1977).  The  equinox  maximum  in  the 
northern  hemisphere  measurements  is  most  pronounced  at  low  lati¬ 
tudes,  with  amplitudes  in  agreement  with  the  southern  hemisphere 
results.  High  latitude  measurements  seem  to  exhibit  a  maximum  in 
winter  (Fukuyama,  1977).  The  emission  rate  also  exhibits  a  small 
nocturnal  variation  of  30%  symmetric  around  a  minimum  at  midnight 
(Kirchoff  et  al.,  198L). 


A  summary  of  the  above  experimental  constraints  is  given  in  Table  2.  We 


will  concentrate  on  the  diurnal  variation  of  the  above  p.nntities. 


seasonal  behavior  which  will  be  useful  for  placing  bounds  on  the  experimental 
values  used  in  comparison  to  our  one-dimens iona 1  model  will  be  discussed  in 


Part  IV. 


4.  DESCRIPTION  OF  MODEL 


The  one-dimensional  model  of  mes  >s;>her  i  c  trace  gases  covers  the  altitude 
range  from  60-120  km.  It  is  an  extension  ot  the  model  previously  used  in 
studies  of  stratospheric  chemistry  (Sze,  Ko,  Specht,  and  Livshits,  1980). 
Concentrations  of  short-lived  photochemical  species  are  calculated  using  the 
grouping  technique  previously  employed  in  stratospheric  models  (Sze  et  al., 
1980). 
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**  '  *  *  ******  *  -  •  -*  -*  J*. 


6 


The  species  calculated  in  our  model  include: 
oxygen  species,  0x:  O^D),  0,  O3 

odd  hydrogen  species,  H0x:  H,  OH,  HO2 ,  H2O2 

odd  nitrogen  species,  N0x:  N,  NO,  NO2,  HNO3,  NO3,  N2O5,  HO2NO2 
sodium  species,  NaX:  Na,  NaO,  Na02 ,  NaOH 
In  addition,  02(1Ag)  is  also  included.  Although  the  N0X  species  and  02(!Ag) 
are  of  minor  importance  in  the  neutral  chemistry  of  the  odd-oxygen  and  odd- 
hydrogen  species,  they  are  included  in  the  chemistry  scheme  in  anticipation  of 
their  roles  in  the  ion  chemistry.  The  effect  of  the  odd-chlorine  and  methyl 
families  on  the  above  species  is  negligible  at  these  altitudes;  thus,  they 
have  not  been  included  in  the  model. 

Background  species  which  are  held  fixed  in  the  model  include  H2,  H2O,  CO, 
N2O ,  O2 ,  and  N2 .  Densities  of  the  major  constituents  (N2  and  O2)  and  tempera¬ 
tures  are  taken  from  the  U.S.  Standard  Atmosphere  (1976).  The  densities  of 
the  background  trace  species  CO,  H2 ,  and  H2O,  are  taken  from  Allen,  Lunine, 
and  Yung  (1984). 

Photochemical  lifetimes  for  H0X  and  0X  are  shorter  than  one  day  at  alti¬ 
tudes  below  80  km.  The  concentrations  of  all  constituents  are  calculated  from 
the  time-dependent  equation 


dn 

dt 


i 


pi  '  LiV 


(3) 


where  n^,  P^,  and  denote  the  density,  production,  and  loss  frequency  for 
the  ic^  species.  Equation  (3)  is  solved  subject  to  periodic  boundary  condi¬ 
tions  over  a  24-hour  period. 


n^  (t  +  24  hrs.)  =  n^(t), 


(4) 


Note  that  the  effect  of  transport  is  ignored  in  this  region. 

Time  constants  for  H0x,  0  ,  and  N0X  become  much  longer  than  one  day  above 
90  km.  The  densities  of  the  above  constituents  are  then  determined  by  down- 
gradient  diffusion.  Although  it  is  within  the  capability  of  the  model  to  cal¬ 
culate  the  distribution  of  the  long-lived  species,  we  have  fixed  their  densi¬ 
ties  in  our  calculations  to  the  profiles  given  by  Allen  et  al.,  (1984)  for  H0x 
and  0  ,  and  by  Solomon,  Crutzen,  and  Roble  (1982)  for  NO  .  Slow  species  are 

A  X 


calculated  In  a  self-consistent  manner  when  we  consider  seasonal  variations  of 
the  vertical  eddy  diffusion  coefficient  (see  Part  IV). 

Sharp  vertical  gradients  in  the  densities  of  0x  and  H0x  result  in  very 

short  (~  1  day)  transport  time  constants  between  80  and  85  km.  The  diurnal 

variation  of  the  above  species  is  then  determined  by  both  photochemical  and 

transport  processes.  The  effects  of  transport  are  approximately  accounted  for 

* 

in  the  present  calculation  by  introducing  an  extra  time-independent  term,  P. 
in  equation  (3).  The  magnitude  and  sign  of  P^  is  chosen  such  that  the  solu¬ 
tions  of  equation  (3)  at  noon  agree  with  the  values  calculated  by  Alien  et 
al.,  (1984)  at  the  same  time. 

The  total  density  of  sodium  species  (NaX)  is  calculated  assuming  an  in¬ 
coming  flux  of  2  x  104  cm  2  s  1  at  100  km.  The  sodium  is  transported  by  dif¬ 
fusion  into  the  troposphere  where  it  is  removed  by  heterogeneous  processes 
with  an  effective  lifetime  of  10  days.  The  values  of  the  diffusion  coeffi¬ 
cients  below  70  km  are  the  same  as  those  used  in  the  stratospheric  model  (Sze, 
Ko,  Specht,  and  Livshits,  1980).  Above  70  km,  a  value  of  106  cm2  s  1  is 
adopted  as  suggested  by  Hunten  (1975).  The  calculated  total  NaX  volume  mixing 
ratio  ranges  from  2.3  x  10  9  at  100  km  to  1  x  10  11  at  60  km. 

The  chemical  scheme  for  the  oxygen,  hydrogen,  and  nitrogen  species  is 
similar  to  the  one  used  in  the  stratospheric  model.  Reaction  rates  are  taken 
from  the  latest  NASA  review  (JPL,  1982),  with  the  corrections  suggested  by 
Allen  et  al.,  (1984).  The  reactions  involving  sodium  species  are  taken  from 
Sze,  Ko ,  Murad,  and  Swider  (1982)  as  given  in  Table  1. 


RESULTS 


Diurnal  profiles  for  the  main  0  species  (0,  O3)  at  82,  84,  86,  and  88  km 


are  shown  in  Figures  2a,  2b,  2c,  and  2d,  respectively;  similar  profiles  at  the 


same  altitudes  are  shown  for  the  HO  species,  H,  OH,  and  HO?  in  Figures  3a, 


3b,  3c,  and  3d.  The  calculations  are  performed  for  38°N,  summer  solstice  con¬ 
ditions.  Our  results  are  similar  to  the  diurnal  variation  calculated  by  Allen 
et  al.,  (1984).  The  Ox  and  H0x  families  are  dominated  by  0  and  H,  respective¬ 
ly.  The  change  in  the  diurnal  behavior  of  these  species  with  altitude  illus¬ 
trate  the  transition  between  the  photocheml cally-control led  region  below  82  km 
and  the  diffusion-controlled  region  above  90  km.  Since  the  sodium  chemistry 
above  00  km  has  very  little  effect  on  the  background  neutral  chemistry,  the 


above  profiles  are  used  for  all  the  different  models  of  sodium  considered  be¬ 
low. 

The  results  shown  in  Figure  2  indicate  that  the  [ 0 ] / [ 0 3 )  ratio  is  102  - 
1 0 3  between  80  and  90  km.  In  the  absence  of  temporary  reservoirs,  [Na]  and 
[NaO]  would  be  controlled  by  the  Chapman  scheme.  The  rates  for  (RRl)  and 
(RR2)  given  in  Table  1  would  then  imply  a  [Na]/[NaO]  ratio  of  101  -  102, 
placing  the  atomic  sodium  peak  well  below  80  km.  Thus,  the  observed  position 
of  the  Na  peak  and  the  constraints  on  (RRl)  and  (RR2)  as  suggested  by  thermal 
chemical  data  strongly  suggest  the  existence  of  temporary  reservoirs  for 
atomic  sodium. 

In  our  previous  study  (Rodriguez  et  al.,  1984),  we  focused  on  the 
interaction  of  the  Na-NaO  species  with  the  temporary  reservoirs,  Na02  and 
NaOH.  In  the  chemical  scheme  illustrated  in  Figure  1,  the  interaction  between 
Na  and  N'a02  is  controlled  by  the  rapid  formation  of  Na02,  (reaction  (RR3)), 
and  recycling  of  Na02  to  Na-NaO  through  (RR6).  Similarly,  coupling  to  NaOH  is 
controlled  by  (RR7)  and  the  photolysis  reaction  (RR4 ) .  At  that  time, 
experimental  data  existed  only  on  (RR3)  (Silver  et  al.,  1984b)  and  (RR4) 
(Rowland  and  Makide,  1982).  Since  then,  (RR7)  has  been  measured  to  have  a 
value  of  2  ±  1  x  10  10  cm3  sec  1  (Ager  and  Howard,  1985).  Given  the 
sensitivity  of  the  model  results  to  (RR7),  we  examine  here  the  sensitivity  of 
results  to  (RR7)  within  the  bounds  of  the  experimental  values  in  conjunction 
with  (RR6). 

Three  models  were  considered,  representing  different  strengths  of  coup¬ 
ling  between  the  Na-NaO  species  and  the  Na02  and  NaOH  temporary  reservoirs. 

The  adopted  parameters  for  these  models  are  listed  in  Table  3.  Model  A  repre¬ 
sents  the  conditions  for  weak  coupling  to  temporary  reservoirs  consistent  with 
measured  kinetic  rates.  The  value  adopted  for  reaction  (RRb)  is  small,  but 
sufficient  to  recycle  enough  Na02  back  to  Na-NaO  so  that  Na  does  not  disappear 
completely  (Sze  et  al.,  1982)  due  to  the  fast  rates  measured  for  reaction 
(RR3).  Model  B  adopts  the  mean  measured  rate  for  the  formation  of  NaOH  (-RR7), 
together  with  small  coupling  to  the  Na02  reservoir.  Model  C  represents 
conditions  of  maximum  coupling  to  the  Na02  reservoir,  with  the  rate  of  (RRb) 
near  the  kinetic  limit. 

Noon  altitude  profiles  for  Na  are  shown  In  Figure  4  for  the  different 
models.  Measurements  by  Siraonich  et  al.  (1979)  for  midnight,  and  Granier  and 


Megie  (1982)  for  noon  are  also  shown  at  altitudes  above  85  km.  As  indicated 
by  the  data,  the  diurnal  variation  of  Na  is  negligible  between  85  and  95  kra. 
Agreement  between  these  measurements  and  Model  C  is  poor.  Better  agreement  is 
obtained  with  Models  A  and  B,  especially  in  the  region  of  the  peak.  The 
discrepancy  between  observation  and  the  calculated  results  above  95  km  is 
probably  due  to  the  effect  of  ion  chemistry.  We  will  discuss  this  in  more 
detail  in  the  Part  Ill. 

The  diurnal  variation  of  the  total  Na  abundance  is  presented  in  Figure  5, 
again  for  all  three  models.  The  absolute  magnitude  of  the  abundance  is  in 
good  agreement  with  observed  values  for  Models  A,  B,  and  C.  The  diurnal 
variation  of  the  total  abundance  is  consistent  with  observations  for  Model  A 
(14%)  but  is  too  large  for  Models  B  and  C  (23%  and  31%,  respectively).  The 
larger  diurnal  variation  in  Model  B  is  a  result  of  increased  abundances  of 
NaOH  due  to  the  increase  in  the  rate  of  (RR7)  by  a  factor  of  2  over  Model  A. 
The  large  diurnal  variations  calculated  for  Model  C  follow  from  the  results 
that  much  of  the  atomic  sodium  is  found  below  82  km  where  the  diurnal 
deviations  of  Na  is  the  largest. 

The  results  presented  in  Figures  4  and  5  suggest  that  Model  C  does  not 
reproduce  the  observed  Na  altitude  profile  and  diurnal  variation  of  the  sodium 
abundance.  Models  A  and  B  yield  very  similar  results  for  the  altitude 
profile.  Based  on  the  diurnal  variation  of  the  column  density,  Model  A  is 
preferable  to  Model  B.  The  diurnal  variation  in  Model  B  is,  however,  close 
enough  to  the  observational  upper  limit,  and  we  cannot  discriminate  between 
values  of  I  and  2  x  10  10  cm3  sec  1  for  the  rate  of  (RR7)  without  considering 
other  second-order  effects  such  as  the  effect  of  uncertainties  in  other 
unknown  rates  or  background  species,  0  and  H2O,  in  particular. 

Figures  6a  and  6b  illustrate  the  diurnal  variation  of  Na  at  82,  84,  and 
86  km,  for  Models  A  (Figure  6a)  and  B  (Figure  6b).  Although  the  diurnal  vari¬ 
ation  in  Model  B  is  slightly  larger  than  in  A,  both  models  yield  results  in 
satisfactory  agreement  with  experimental  data  below  85  km.  We  should  note 
that  the  day-to-night  ratio  changes  dramatically  over  an  altitude  range  of  2 
km  in  this  region.  This  fact  should  be  kept  in  mind  in  comparing  theoretical 
results  to  lidar  observations,  since  the  latter  typically  have  an  altitude 
resolution  of  1-2  km  (Megie  and  ULamont,  1977;  Simonich  et  al.,  1979). 


The  interaction  between  the  different  sodium  reservoirs  is  illustrated  in 
Figures  7a  through  7d.  In  the  region  80-90  km,  the  [ Na  J / [ Na02 ]  ratio  can  be 
estimated  from  the  approximate  chemical  equilibrium  conditions  for  Na02: 

[Na]  [02  ]  [M]  =  k6[Na02][0].  (5) 

In  both  Models  A  and  B,  Na  is  the  most  abundant  species  at  90  km  during  the  day, 
with  NaOH  comparable  to  Na  at  night  for  Model  B.  At  lower  altitudes,  with  larger 
values  for  [O2 ]  and  [M]  and  smaller  values  for  [0],  Na02  becomes  comparable  to 
NaOH  at  82  km.  The  [ NaOH ] / [ Na02 ]  ratio  can  be  obtained  from  the  approximate 
re lat ions 

(J4  +  k12  [H])  [NaOHj  =  f  Na02  ]  (k1Q[H]  +  kn[0H])  +  k7[Na0][H20]  (6) 

kj [Na][03 ]  =  k2 [ NaO ] [ 0]  (7) 

and  equation  (5).  Rapid  changes  in  [NaOH]  are  expected  at  dawn  and  dusk  due  to 
the  onset  and  disappearance  of  photolysis. 

The  diurnal  variation  of  Na  can  be  related  to  the  variation  in  Na02.  0,  and 
NaOH  through  (5)  -  (7).  At  90  km,  atomic  oxygen  exhibits  no  diurnal  variation, 
and  changes  in  Na02  and  Na  are  due  to  rapid  changes  in  NaOH  from  day  to  night. 

As  we  increase  the  rate  of  (RR7)  from  Model  A  to  Model  B,  the  Na()H/Na02  ratio 
increases,  and  the  large  diurnal  changes  in  NaOH  induce  correspondingly  larger 
variations  in  Na02  and  Na  near  the  sodium  layer  peak  at  90  km.  On  the  other 
hand,  atomic  oxygen  exhibits  a  large  diurnal  variation  at  82  km  and  determines 
the  diurnal  variation  of  Na  to  first  order  at  that  altitude. 

b.  DISCUSSION  AND  CONCLUSIONS 

Results  of  Models  A,  B,  and  C  are  summarized  in  Table  4.  The  experimental 
results  from  Table  2  are  included  for  comparison.  Also  shown  in  Table  4  is  the 
midnight  D-line  emission  rate.  We  note  that  the  calculated  D-iine  emission  rate 
for  a  given  model  depends  on  the  assumed  branching  ratio  for  production  of  Na(2P) 
by  reaction  (RR2),  for  which  only  estimates  exist  (Sze  et  al.,  1982;  Bates  and 
Ujha,  1980).  We  have  adopted  a  value  of  .3  for  the  branch  producing  Na(2P). 
Results  on  the  D-line  emission  rate  are  reasonable  for  all  models,  and  thus  are 
inconclusive  for  our  purposes.  As  noted  previously.  Model  C  agrees  poorly  with 
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observations  of  the  peak  height,  peak  density,  the  column  density,  and  diurnal 
variation.  We  can  therefore  ruLe  out  Model  C  .is  an  unrealistic  description  ot 
sodium  chemistry  in  the  upper  atomsphere.  The  rate  for  (RRb)  (Na02  +  0  NaO  + 

C>2  )  must  therefore  be  significantly  lower  than  10~10  cm3  s  3,  closer  to  10  13  cm3 

s'1  . 

It  is  harder  to  draw  any  definite  conclusions  for  differentiating  between 
Models  A  and  B.  The  best  basis  for  differentiation  is  the  diurnal  variation  of 
local  densities  and  total  column  densities  of  Na.  Local  density  variations  for 
both  models  are  consistent  with  observations,  given  the  altitude  resolution  of 
lidar  measurements.  Variations  in  the  total  column  density  for  Models  A  and  B 
are  more  sharply  differentiated,  and  the  question  thus  lies  on  the  interpretation 
of  the  experimental  data. 

Early  lidar  measurements  of  the  sodium  layer  yielded  day-to-night  ratios  of 
at  most  10%  (Gibson  and  Sandford,  1972)  or  changes  during  the  night  of  about  4% 
(Simonich  et  al.,  1979).  More  recent  observations  of  Clemesha  et  al.  (1982) 
which  were  carried  out  on  a  continuous  basis  yielded  diurnal  variations  with  an 
amplitude  of  ~  13%  from  the  mean  (30%  peak-to-peak) .  The  observations,  however, 
exhibited  a  strong  semidiurnal  component,  and  the  above  authors  suggested  that 
the  diurnal  variation  is  due  entirely  to  tidal  effects.  Thus,  definite 
conclusions  could  ne  made  only  after  incorporation  of  tidal  effects  into  existing 
photochemical  models,  coupled  with  higher  resolution  measurements  of  the  diurnal 
variation  of  Na.  We  thus  cannot  constrain  the  rate  (RR7)  (NaO  +  H20  +  NaOH  +  OH) 
to  more  than  a  factor  of  two  at  this  point.  The  lower  bound  used  in  Model  A  is 
consistent  with  measurements  by  Ager  and  Howard  (1985). 

Another  consequence  of  the  high  rate  for  reaction  (RR3)  is  that  the  diurnal 
variation  of  atomic  sodium  below  the  peak  can  be  estimated  to  first  order  in 
Models  A  and  B  from  the  photochemical  equilibrium  relation  between  Na  and  Na02 
given  in  equation  (5).  The  diurnal  variation  of  Na  then  follows  that  of  atomic 
oxygen  if  NaOp  does  not  change  appreciably,  a  condition  satisfied  below  the  Na 
peak  in  Model  A .  Increased  coupling  to  NaOH  leads  to  a  higher  diurnal  variation 
ot  Na02  in  Model  B.  Comparison  of  the  diurnal  variation  of  0  in  Figure  2  and 
that  ot  Na  in  Figure  f>  suggests  general  agreement  between  both  diurnal  changes. 
The  ratio  I ()  1  pus  /  i 0 1  |)lwn  r>*3  at  82  km  anf*  1-7  at  84  km.  The  additional 
coupling  to  Na‘>i|  is  reflected  in  larger  ratios  for  all  models,  particularly  at 
84  km.  In  Model  C,  the  diurnal  variation  ot  atomic  oxygen  is  no  longer  me 


first-order  process  contributing  to  the  variability  of  sodium,  even  at  82  km.  In 
this  case,  the  abundance  of  Na02  is  much  less  than  that  of  NaOH,  and  the  diurnal 
variation  of  Na  is  determined  by  the  interaction  with  the  NaOH  reservoir.  As 
stated  above,  however,  we  feel  that  this  model  does  not  represent  the  chemistry 
of  sodium  species. 

Given  that  the  diurnal  variation  of  atomic  sodium  depends  crucially  on  the 
diurnal  variation  of  atomic  oxygen,  the  observed  absence  of  appreciable  diurnal 
change  in  the  total  Na  abundance  requires  that  the  bulk  of  the  sodium  be  located 
in  a  region  where  0  is  constant  throughout  the  day,  i.e.,  about  85  km.  Once  the 
total  NaX  is  fixed,  the  position  of  the  peak  is  determined  to  first  order  by  the 
magnitude  of  (RR3).  Lowering  the  peak  by  modifying  the  NaX  profile  through 
changes  in  the  Na  source  or  eddy  diffusion  coefficient  could  be  reflected  in  a 
larger  diurnal  variation  of  the  column  density  of  sodium.  Since  the  time 
constant  for  the  Na~Na02  exchange  ranges  from  seconds  to  a  few  minutes  below  90 
km,  any  modelling  of  transport  effects  (i.e.,  tides,  gravity  waves)  must  include 
a  self-consistent  treatment  of  the  atomic  oxygen  variation. 

The  photolytic  reactions, 

NaOH  +  hv  -►  Na  +  OH,  (RR4) 

Na02  +  hv  >  Na  +  O2 ,  (RR5) 

could  also  play  an  important  role  in  the  diurnal  variation  of  sodium  if  their 
values  were  very  different  from  those  assumed  here.  Reaction  (RR5)  could  com¬ 
pete  with  the  reaction  with  0  (RR6)  if  J5  ~  10  2  s  l.  In  such  a  case,  the 
recycling  of  Na02  to  Na  would  occur  via  reaction  (RR5),  and  the  constraints 
derived  here  for  reaction  (RR6)  can  be  easily  translated  into  similar 
constraints  for  reaction  (RR5).  Photolysis  of  Na02  is  also  important  in  the 
stratosphere,  where  the  concent ra t Ion  of  atomic  oxygen  is  small. 

The  value  adopted  for  reaction  (RR4)  by  Rowland  and  Makide  (1982)  assumes 
that  all  the  observed  EUV  absorption  of  NaOH  is  due  to  photolysis.  If  this  is 
not  the  case,  and  J*  is  smaller,  our  models  would  produce  more  NaOH,  but  its 
day-to-night  change  would  be  smaller  and  the  effects  of  coupling  of  Na-NaO  to 
the  NaOH  reservoir  would  be  consequently  reduced. 

Ionization  of  atomic  sodium  can  also  affect  our  results  in  the  region 
above  the  peak.  Atomic  sodium  is  ionized  both  by  photons  and  by  change  ex¬ 
change  with  02+,  with  an  effective  rate  of  ~  10  5  s  1  (Swider,  1970;  1984). 

Loss  of  Na+  occurs  primarily  through  clustering  with  N?, 
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N2  +  M 


(7) 


Na+  +  N2  +  M  ->•  Na+  • 

foLlowed  by  dissociative  recombination  of  Na+  •  N2  (Richter  and  Sechrist, 
1979).  The  partitioning  between  Na  and  Na+  could  thus  be  determined  by  the 
rate  of  the  clustering  reaction  (7),  which  has  not  been  measured.  Vertical 
ion  drifts  could  also  change  the  Na/Na+  ratio.  This  would  reduce  the  amount 
of  atomic  Na  calculated  in  the  model,  bringing  the  result  into  much  better 
agreement  with  the  measurements  shown  in  Figure  4  above  the  peak.  We  will 
discuss  this  matter  in  more  detail  in  Part  III. 

In  conclusion,  the  results  of  our  study  point  out  the  importance  of  the 
Na02  temporary  reservoir  in  determining  the  diurnal  variation  of  atomic  sodi¬ 
um.  The  diurnal  variation  of  atomic  oxygen  thus  plays  a  dominant  role  in  de¬ 
termining  diurnal  changes  in  Na.  Present  experimental  data  constrains  the 
recycling  of  NaC>2  through  reaction  with  0  (RR6),  limiting  the  rate  to  consi¬ 
derably  less  than  10  13  cm3  s  1 .  Large  amounts  of  NaOH  are  produced  through 
the  reaction  of  NaO  with  H2O  (reaction  (RR7)),  although  the  rate  for  this 
process  has  a  factor  of  3  uncertainty.  Further  clarification  of  this  problem 
requires  measurements  of  unknown  reaction  rates  (in  particular,  (RR4),  (RR5) , 
and  (RR6)),  continuing  high-resolution  observations  of  sodium  species  in  the 
upper  atmosphere,  and  further  modelling  incorporating  tidal  and  gravity  wave 
effects  on  both  sodium  species  and  background  trace  gases. 


Table  l:  Reaction  Scheme  and  Rate  Constants 


Reaction 

Ra  t  e  Coe  1  t i c i e n ta 

1 . 

N'a  +  O3 

-► 

NaO  +  02 

3. 5  x  lo"10 

lb] 

2. 

NaO  +  0 

-► 

Na  (2P,  2S)  +  02 

4  (  —  1  1  ) 

lo] 

3. 

Na  +  On  +  N2 

> 

Na02  +  N2 

1.9  X  1 0_3 0 ( r/300)_1  • 

Md] 

4. 

NaOH  +  hv 

-► 

Na  +  OH 

2  (-3) 

le] 

5. 

Na02  +  hv 

-► 

Na  +  02 

1  ( -4 ) 

[  r  1 

b . 

Na02  +  0 

> 

NaO  +  02 

l(-10) ,  1  ( —  1 3 ) 

Ic] 

7. 

NaO  +  H2  0 

> 

NaOH  +  OH 

(2  ±  1)  x  lo"  10 

[h] 

8. 

NaO  +  O3 

-* 

Na02  +  02 

1 .8  x  10"  10 

[hi 

9. 

NaO  +  H 

-► 

Na  +  OH 

1  (  —  1  4  ) 

Is! 

10. 

Na02  +  OH 

-► 

NaOH  +  02 

1  (  —  1 1 ) 

It  1 

1 1 . 

Na02  +  H 

NaOH  +  0 

K-13) 

IbI 

12. 

NaOH  +  H 

Na  +  H20 

1 . 4  ( - 1 2 ) 

It] 

13. 

NaOH  +  C1 0) 

NaO  +  OH 

1  (-10) 

Igl 

14. 

Na  +  H02 

> 

NaOH  +  0 

1  ( - 1 0 ) 

It] 

15. 

NaO  +  (33 

•+ 

Na  +  202 

1  ( —  1 0 ) 

1  h  | 

1  b . 

NaO  +  H02 

-► 

NaOH  +  02 

K-ii) 

If] 

17. 

NaO  +  H2 

•> 

Na  +  H2  0 

1  ( - 1 1 ) 

Ih] 

18. 

NaO  +  H2 

NaOH  +  H 

2 .  b  (  - 1 1 ) 

Ih] 

Read 

3(-10) 

as  3 

x  10  10. 

Units  are  cm^ 

s  1  , 

except  cmh  ; 

n  1  for  (3)  and 

s"l 

for  (4), 

(5), 

and  (  15). 

] b ] ,  Zahniser  et 

al.,  1985). 

[ c ] ,  Sxe  et 

al . 

(1982). 

[d  ]  , 

Silver  et 

al.  (1984). 

le], 

Rowland  and 

Makide 

(1982).  [f] 

,  Liu 

and  Reid 

(1979).  |gl. 

Kirchoff  and  Cletnesha  (1985). 

Ih] 

Age r  and 

Howard  (1983). 

Table  2:  Experimental  Data  on  Sodium 


References 


Sodium  Layer 

Peak  Altitude 

90  -  92  km 

Megie  and  Blamont,  1977; 

Peak  Density 

3-5  x  103  cm-3 

Simonich  et  al.,  1979 

Total  Abundance 

4  x  109  cm  3 

Gibson  and  Sandford,  1972; 

(2-8  x  109  cm-2) 

Clemesha  et  al.,  1982 

Diurnal  Variation 

Total  Abundance3 

<  15% 

Gibson  and  Sandford,  1971; 
Clemesha  et  al.,  1982 

t Na  1  Dus^  Dawn 

2  (83  km) 

Simonich  et  al.,  1979; 

10  (81  km) 

Clemesha  et  al.,  1982 

D-line  Emission  Rate 

~  50  R 

Fukuyama,  1977;  Kirchoff  et 

(30-100  R) 

al.,  1981 

Percent  deviation  from  median  abundance 


Table  3:  Adopted  Models 


k  (NaO?  +  0) 


( cm3  s  1  ) 


k  (NaO  +  H?0) 


(cm 3  s  1 ) 


Model  A 
Model  B 
Model  C 


2  x  10  10 


Table  4:  Calculated  Sodium  Parameters  for  the  Different  Models 


Model  B 


Model  A 


Model  C 


Measured  Value 


Peak  Altitude 


Peak  Density 


90  km 


90  km 


86  km 


90-92  km 


4.3  x  103  cm  3  4.2  x  103  cm  3  6.3  x  103  cm  3  3-5  x  103  cm  3 


Total  Abundance3  4.47  x  10<3  cm  2  4.13  x  103  cm  2  6.34  x  103  cm  2  2-8  x  103  cm  2 


Diurnal  Variation, 
Total  Abundance 


^Wsk^^Dawn 


D-Line  Emission 
Rate  (night) 


14% 


7.1  (82  km) 

2.4  (84  km) 


35  R 


23% 


7.4  (82  km) 
3.1  (84  km) 


46  R 


31% 


25.5  (82  km) 
3.8  (84  km) 


47  R 


<15;, 


10  (81  km) 
2  (83  km) 


50  R 


^Averaged  over  one  day. 

3Percent  deviation  from  the  median  abundance. 
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Figure  1.  Chemical  reactions  controlling  the  behavior  of  sodium  species  in 
the  mesosphere.  The  dotted  box  in  the  left-hand  side  represents 
the  scheme  proposed  by  Chapman  (1939).  Coupling  to  the  temporary 
reservoirs  NaC>2  and  NaOH  occurs  primarily  through  the  reactions 
denoted  by  solid  arrows.  The  dotted  arrows  represent  direct 
production  of  NaOH  from  Na02  (Liu  and  Reid,  1979).  This  mechanism 
is  the  primary  source  of  NaOH  if  the  reaction  rate  of  NaO  with  H 2O 
is  negligible. 

Figure  2.  Diurnal  variation  of  odd-oxygen  species  0  and  O3,  calculated  by  our 
model  at  (a)  82  km,  (b)  84  km,  (c)  86  km,  and  (d)  88  km.  Input 
parameters  and  assumptions  of  our  model  are  discussed  in  the  text. 

Figure  3.  Diurnal  variation  of  odd-hydrogen  species  H,  OH,  and  HOX  calculated 
at  (a)  82  km,  (b)  84  km,  (c)  86  km,  and  (d)  88  km. 

Figure  4.  Altitude  profiles  of  atomic  sodium  at  noon  calculated  for  Models  A, 
B,  and  C.  The  measurements  of  Simonich  et  al.  (1979)  at  midnight 
(x)  and  of  Granier  and  Megie  (1982)  at  noon  (o)  are  also  included 
for  comparison. 

Figure  5.  Diurnal  variation  of  the  total  column  abundance  of  atomic  sodium, 
calculated  for  Models  A,  B,  and  C. 

Figure  6.  Diurnal  variation  of  local  densities  of  atomic  sodium  at  82,  84, 
and  86  km.  Results  are  shown  for  Model  A  (6a)  and  Model  B  (6b). 

Figure  7.  Diurnal  variation  of  the  local  densities  of  sodium  species  at  82 
and  90  km,  calculated  for  Model  A  (7a  and  7b)  and  Model  B  (7c  and 
7d). 
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Chemical  reactions  controlling  the  behavior  of  sodium  species  in 
the  mesosphere.  The  dotted  box  in  the  left-hand  side  represents 
the  scheme  proposed  by  Chapman  (1939).  Coupling  to  the  temporary 
reservoirs  Na02  and  NaOH  occurs  primarily  through  the  reactions 
denoted  by  solid  arrows.  The  dotted  arrows  represent  direct 
production  of  NaOH  from  Na02  (Liu  and  Reid,  1979).  This  mechanism 
is  the  primary  source  of  NaOH  if  the  reaction  rate  of  NaO  with  H2O 
is  negligible. 
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Figure  2.  Diurnal  variation  of  odd-oxygen  species  0  and  O3  ,  calculated  by  our 
model  at  (a)  82  ka,  (b)  84  km,  (c)  86  kra,  and  (d)  88  km.  Input 
parameters  and  assumptions  of  our  model  are  discussed  In  the  text. 
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Altitude  profiles  of  atomic  sodium  at  noon  calculated  for  Models  A 
B,  and  C.  The  measurements  of  Simonich  et  al.  (1979)  at  midnight 
(x)  and  of  Granier  and  Megle  (1982)  at  noon  (o)  are  also  included 
for  comparison. 
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PART  II 

Possible  Iapact  of  Sodium  Species  on  CIO 
and  O3  in  the  Upper  Stratosphere 

ABSTRACT 


Recent  measurements  of  the  rates  for  reactions  of  Na02  and  NaOH  with  HC1 
and  of  the  photoabsorption  cross  section  of  NaCl  allow  for  realistic  estimates 
of  the  impact  of  sodium  chemistry  on  chlorine  partitioning  in  the  upper  strat¬ 
osphere.  When  sodium  chemistry  is  incorporated  into  our  one-dimensional  model 
of  the  upper  atmosphere,  the  results  indicate  that  the  calculated  CIO  concen¬ 
trations  near  50  km  could  increase  by  up  to  a  factor  of  two  with  a  correspon¬ 
ding  decrease  in  calculated  O3  abundances  of  as  much  as  15  percent.  The  im¬ 
pact  on  CIO  and  O3  depends  crucially  on  the  assumed  meteoritic  source  of  sodi¬ 
um,  on  transport  from  the  mesosphere,  and  on  transformation  and  removal  of 
gaseous  sodium  in  the  stratosphere.  Since  the  fluctuations  in  the  abundance 
of  gaseous  sodium  could  induce  O3  variations  comparable  to  those  caused  by 
changes  in  temperature  and  solar  UV,  consideration  of  the  chemistry  of  sodium 
and  other  metals  may  be  an  important  element  in  efforts  toward  early  detection 
of  trends  in  upper  stratospheric  O3 .  Future  increases  in  stratospheric  chlor¬ 
ine  will  magnify  the  impact  of  sodium  species  on  ozone. 
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i  . 

The  catalytic  efficiency  of  chlorine  compounds  in  removing  stratospheric 
O3  is  particularly  sensitive  to  the  abundance  of  active-chlorine  species  (Cl, 
CIO).  The  chlorine  catalytic  cycle  is  most  effective  in  the  present  atmo¬ 
sphere  between  40  and  50  km,  representing  about  one  third  of  the  total  loss  of 
O3  at  these  altitudes.  It  is  thus  in  this  region  where  an  early  detection  of 
ozone  perturbations  due  to  anthropogenic  chlorine  is  most  probable. 

Our  present  understanding  of  chlorine  chemistry  indicates  that  the  parti¬ 
tion  between  active-chlorine  and  the  reservoir  species  HC1  is  primarily  deter- 


mined  by  the  reactions 

OH  +  HCl  ->■  Cl  +  H20 

(1) 

Cl  +  CH4  >  HCl  +  CH3 

(2) 

Cl  +  H02  +  HCl  +  02 

(3) 

Other  processes  which  could  affect  the  above  partition  include  reactions  of 
the  type 

MOH  +  HCl  >  MCI  +  U20  (4) 

M02  +  HCl  +  MCI  +  HO2  (5) 

where  M  denotes  stratospheric  metal  species  of  raeteoritic  origin  (Na,  Ca,  Ni, 
Mg,  Al,  Si,  and  Fe,  among  others).  The  potential  importance  of  reaction  (4) 
was  first  recognized  by  Murad  et  al.  (1981)  who  suggested  that  the  species  MCI 
could  act  as  an  important  reservoir  of  chlorine  if  its  stratospheric  abundance 
is  comparable  to  that  of  HCl.  On  the  other  hand,  reactions  (4)  and  (5)  fol¬ 
lowed  by  photolysis  of  MCI 

MCI  +  hv  >  M  +  Cl  (6) 

could  constitute  a  catalytic  cycle  for  recycling  HCl  to  active  chlorine  (Row¬ 
land  and  Rogers,  1481). 

ReaListic  estimates  of  the  effect  of  the  above  reactions  have  been  ham¬ 
pered  by  a  lack  of  kinetic  data  for  metal  species.  Recent  measurements  of  the 
rates  of  reactions  of  NaOH  and  NaO?  with  HCl  indicate  very  fast  rates  ap- 


proaching  the  kinetic  limit  (Silver  et  al.,  1984a;  Silver  and  Kolb,  1986). 
These  kinetic  data  can  be  combined  with  the  latest  measurement  of  the  photoly¬ 
sis  rate  of  NaCl  (Silver  et  al.,  1986)  to  provide  a  realistic  assessment  of 
the  importance  of  (4)  to  (6)  in  the  recycling  of  chlorine  species  by  sodium. 

We  report  here  results  of  such  calculations. 

2.  CHEMISTRY  OF  SODIUM  SPECIES 

Direct  observations  of  atomic  sodium  in  the  Earth's  mesosphere  have 
established  the  existence  of  a  narrow  layer  of  Na  with  average  peak  densities 
of  2-5  x  1()3  cm  3  near  90  km  (cf.  Clemesha  et  al.,  1982).  No  other  sodium 
compounds  have  been  observed  directly.  Their  atmospheric  abundance, 
particularly  in  the  stratosphere,  must  therefore  be  inferred  from  models  in¬ 
corporating  sources,  chemistry,  and  transport  of  Na  species  and  constrained  by 
the  mesospheric  observations. 

We  have  incorporated  the  chemistry  of  active  sodium  species  NaX  (Na,  NaO, 
Na02 ,  NaOH,  NaCl)  into  the  A£R  1-D  model  (Ko  and  Sze,  1984),  now  covering  the 
altitude  range  0-120  km.  Adopted  kinetic  and  photochemical  data  are  listed  in 
Table  l,  and  the  chemical  scheme  is  illustrated  in  Figure  1.  Reactions  R1  and 
R2  modulate  rapid  conversion  between  Na  and  NaO  near  90  km  and  account  for  the 
sodium  D-line  nightglow  (Chapman,  1939).  Reaction  of  Na  with  O2  and  NaO  with 
H2O  (R3  and  R7)  sequester  most  of  the  Na  and  NaO  in  the  temporary  reservoirs 
Na02  and  NaOH  between  60  and  85  km.  The  above  reservoir  species  are  in  turn 
converted  to  NaCl  below  60  km  through  reaction  with  HC1.  The  impact  of 
chlorine  chemistry  on  the  sodium  layer  at  90  km  is  negligible. 

Sodium  is  assumed  to  be  produced  in  atomic  form  by  meteor  ablation  be¬ 
tween  80  and  90  km  (Hunten,  1981).  Estimates  of  the  average  influx  range  from 
5  x  10^  to  4  x  104  cm  2  s  1  (Kirchoff  and  Clemesha,  1983;  Hunten,  1981) 
depending  primarily  on  the  adopted  mesospheric  eddy  diffusion  coefficient  and 
Na  density  at  90  km.  We  adopt  in  our  model  the  shape  of  the  production  pro¬ 
file  of  Hunten  (1981)  and  normalize  it  to  a  total  influx  of  2  x  104  cm  2 
s‘l .  The  primary  loss  for  active  sodium  is  assumed  to  be  washout  in  the 
troposphere  with  a  lifetime  of  ID  days.  Eddy  diffusion  coefficients  in  the 
mesosphere  are  those  of  Allen  et  al.  (1981).  Kinetic  data  and  boundary 
conditions  for  other  stratospheric  trace  gases  are  those  used  in  the  modeling 
of  present-day  conditions  (Watson,  1986). 


Tne  model  described  above  calculates  peak  densities  of  Na  of  4  x  10  3  cm  3 
at  40  km,  and  a  diurnal  variation  of  the  column  density  of  Na  of  less  than  10 
percent,  in  agreement  with  observations  (Clemesha  et  al.,  1982).  Diurnally- 
averaged  profiles  for  the  major  stratospheric  NaX  species  (NaOH,  Na02,  N'aCl) 
calculated  by  this  model  are  shown  in  Figure  2.  Sodium  chloride  is  the  major 
form  of  NaX  below  60  km,  exhibiting  densities  between  10s  and  10°  cm  3.  Such 
abundances  (a  few  pptv)  are  small  in  comparison  to  those  of  HC1  and  CIO  (>  100 
pptv),  implying  a  negligible  role  for  NaCl  as  a  temporary  reservoir  of 
chlorine . 

Figure  1  also  indicates  other  hypothetical  net  sinks  for  NaX  in  the 
stratosphere.  These  sinks,  denoted  by  include:  formation  of  stable  so¬ 

dium  bicarbonate  (Murad  and  Swider,  1979),  attachment  to  aerosol  and  smoke 
particles  (Hunten,  1981),  and  formation  of  stable  NaCl  polymers  (Lamb  and 
Benson,  1986).  The  sensitivity  of  our  calculated  NaX  to  the  polymerization 
process  is  illustrated  in  Figure  5.  The  soLid-line  profile  assumes  tropo¬ 
spheric  washout,  to  be  the  onLy  sink  of  active  sodium,  as  in  Figure  2.  The 
calculations  denoted  by  the  dashed  Line  adopt  the  estimates  of  Lamb  and  Benson 
(198o)  for  the  dimerization  rate  of  NaCl,  as  well  as  their  polymerization 
scheme.  The  dash-dot  line  represents  NaX  profiles  calculated  with  a 
dimerization  rate  a  factor  of  10  smaller,  thus  allowing  for  uncertainties  in 
the  estimate  and  for  possible  reductions  in  the  effective  sink  of  NaX  through 
photolysis  of  NaCl  polymers  or  other  chemical  reactions.  The  relative 
partitioning  of  NaX  species  for  all  the  above  cases  is  the  same  as  in  Figure 
2.  The  calculated  NaX  near  40  km  can  change  by  up  to  two  orders  of  magnitude, 
depending  on  the  adopted  polymerization  mechanism. 

1.  IMPACT  OF  SODIUM  CHEMISTRY  ON  CHLORINE  SPECIES 

The  effectiveness  of  the  catalytic  cycle  (4)  -  (6)  in  converting  HC1  to 
ci  depends  crucially  on  the  abundance  of  NaCl  and  its  photolysis  rate.  Figure 

•  compares  the  calculated  rate-limiting  step  for  this  cycle  with  the 
unvent  >  fi.il  re  eye  1  ing  of  HC  L  through  re...  t  ion  (1).  The  photolysis  rates, 

which  'live  be*  i  calcul  tted  I  ruin  the  UV  cross  sections  recently  measured  by 
.stiver  **t  al.  (1986,),  ire  smaller  than  the  estimates  of  Rowland  and  Rogers 

*  1982).  The  rates  tor  the  above  mechanism  are  comparable  near  60  km. 


The  effect  of  sodium  chemistry  on  the  calculated  CIO  is  shown  in  Figure 
5.  Percent  changes  are  calculated  relative  to  the  CIO  profiles  obtained  with 
no  sodium  chemistry.  The  solid,  dashed,  and  dash-dot  profiles  represent  the 
same  assumptions  for  the  active-sodium  sink  as  in  Figure  3.  We  have  also 
included  in  this  figure  a  case  with  the  NaX  concentrations  a  factor  of  2 
higher  than  in  Figure  2  and  no  polymerization  sink.  Given  the  uncertainties 
in  the  meteoritic  input,  the  mesospheric  vertical  transport,  and  the  observed 
peak  Na  densities,  such  concentrations  of  NaX  cannot  be  ruled  out.  Changes  in 
the  CIO  calculated  above  40  km  range  from  a  few  percent  to  almost  a  factor  of 
two,  depending  on  the  assumed  net  sink  of  NaX. 

The  changes  in  the  calculated  CIO  are  reflected  in  the  calculated  O3 
above  40  km.  Percent  changes  in  the  calculated  O3  are  shown  in  Figure  6,  for 
the  same  cases  as  in  Figure  5.  Maximum  changes  occur  between  50  and  55  km. 

The  calculated  ozone  could  be  reduced  by  as  much  as  10-15  percent  at  these 
altitudes,  but  the  actual  reduction  will  again  be  very  sensitive  to 
assumptions  about  the  polymerization  sink. 

4.  DISCUSSION  AND  CONCLUSIONS 

The  calculations  presented  above  indicate  that  sodium  chemistry  may  play 
an  important  role  in  determining  the  O3  abundance  in  the  upper  stratosphere  by 
influencing  the  partitioning  of  chlorine  species.  Variations  in  the 
stratospheric  NaX  abundances  are  expected  due  to  seasonal  changes  in 
mesospheric  transport  (Garcia  and  Solomon,  1985),  other  dynamical  processes  or 
to  sporadic  changes  in  the  sources  and  sinks.  Sodium  chemistry  could  then 
induce  variations  of  up  to  a  factor  of  two  in  the  CIO  and  up  to  15  percent  in 
the  ozone  at  50  km.  These  changes  in  ozone  are  comparable  to  variations  of 
order  5-10  percent  due  to  dynamical,  temperature,  or  solar  cycle  variations 
(cf.  Keating  et  al.,  1986)  and  would  be  even  larger  for  an  enhanced 
stratospheric  chlorine  burden.  The  possible  effect  of  sodium  chemistry  on 
ozone  variations  should  thus  be  included  in  efforts  towards  early  detection  of 
the  perturbation  to  upper  stratospheric  ozone  by  increases  in  chlorine. 

The  calculated  peak  density  of  the  sodium  layer  at  90  km,  which  is  the 
only  observational  constraint  on  stratospheric  Na  abundances,  depends  not  only 
on  the  adopted  Na  influx  but  also  on  the  assumed  eddy  coefficient  near  the 
sodium  peak  (Liu  and  Reid,  1979)  and  on  the  relative  magnitude  of  the  rates  of 
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reaction  of  Na  with  O2  (R3)  and  of  Na(>2  with  0  (R6)  (Rodriguez  et  al., 

1984).  Better  estimates  of  the  stratospheric  sodium  abundance  require  better 
determinations  of  the  flux,  composition  and  ablation  of  meteorites  and 
analysis  of  observations  of  mesospheric  Na,  O3,  0,  and  sodium  D  line  emission 
in  conjunction  with  improved  kinetic  data  on  sodium  and  realistic  models  of 
the  chemistry  and  dynamics  of  the  mesosphere.  Measurements  of  the 
polymerization  rate  of  NaCl  and  of  homogeneous  and  heterogeneous  processes  af¬ 
fecting  gaseous  sodium  are  also  needed. 

Observational  validation  of  the  calculations  presented  above  could  best 
be  obtained  by  detection  of  gaseous  NaCl  in  the  upper  stratosphere. 
Observations  of  O3  and  CIO  above  50  km  could  also  be  analyzed  to  detect 
possible  seasonal  changes  induced  by  seasonal  variations  in  the  transport  of 
Na  from  the  mesosphere  by  the  mechanism  of  Garcia  and  Solomon  (1985). 

Our  calculations  also  suggest  the  possible  importance  of  other  metals  of 
meteoritic  origin.  The  raeteoritic  abundance  of  metals  such  as  Ca  and  Al  are 
comparable  to  that  of  Na,  while  those  of  Si,  Mg  and  Fe  could  be  up  to  one 
order  of  magnitude  higher  (Mason,  1971).  Very  little  is  known  about  the  chem¬ 
istry  and  atmospheric  distribution  of  the  above  species,  although  we  may  ex¬ 
pect  similar  chemical  schemes  as  that  of  sodium  (Murad  and  Swider,  1979).  The 
results  presented  here  thus  warrant  further  investigation  of  the  role  of  these 
metals  in  the  chemistry  of  the  upper  stratosphere. 
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Table  1 

Reaction  Scheme  and  Rate  Constants 


Reaction 

Rate3 

Rl. 

Na  +  O3  -*■  NaO  +  O2 

3. 5(-10)b 

• 

CM 

at 

NaO  +  0  -*•  Na  (2P,  2S)  +  02 

4 ( — 1 1 )h* 

R3. 

Na  +  O2  +  N2  ■*  Na02  +  N2 

x  (T/300)”1 ’ ld 

1.9  x  10_30 

R4. 

NaOH  +  HC1  -*•  NaCl  +  H20 

2.0(-10)c 

R5. 

Na02  +  HC1  -*•  NaCl  +  H02 

2.3(-10)e 

R6. 

Na02  +  0  +  NaO  +  02 

1 ( — 1 3 ) h* 

R7 . 

NaO  +  H20  NaOH  +  OH 

2(-10)g 

R8. 

NaO  +  O3  Na02  +  02 

1.8(-10)g 

R9 . 

Na02  +  OH  -*•  NaOH  +  02 

l(-ll)h* 

RIO. 

Na02  +  H  +  NaOH  +  0 

1 ( —14 )h* 

Rll. 

NaO  +  O3  +  Na  +  202 

l(-10)g 

R12 . 

NaO  +  H2  -*•  Na  +  H20 

l(-ll)g 

R13. 

NaO  +  H2  ->•  NaOH  +  H 

2.6(-ll)g 

R14 . 

NaOH  +  hv  -►  Na  +  OH 

2(-3)i 

R15. 

Na02  +  hv  +  Na  +  02 

1 (-4 )h* 

R16. 

NaCl  +  hv  -*■  Na  +  Cl  * 

see  text^ 

aRead  3 ( — 1 0 )  as  3  x  10  10 .  Units  are  cm3  s_1 ,  except  cm6  s  1  for  (3)  and  s  1 
for  (14),  (15),  and  (16).  bZahniser  et  al.  (1985).  cSilver  et  al. 

(1984a).  ^Silver  et  al.  (1984b).  eSilver  and  Kolb  (1986).  ^Silver  et  al. 
(1986).  gAger  and  Howard  (1985).  bSze  et  al.  (1982).  ^Rowland  and  Makide 
(1982).  *Estimated  rates. 
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Chemistry  of  sodium  species  in  the  stratosphere  and  mesosphere. 
Rates  for  processes  denoted  by  dotted  lines  have  not  been  measured. 


Calculated  diurnally-averaged  densities  of  sodium  species  in  the 
upper  stratosphere.  The  species  Na02  and  NaOH  disappear  at  night  in 
this  altitude  range,  while  NaCl  exhibits  little  diurnal  variation. 
The  only  loss  of  NaX  is  by  washout  in  the  troposphere. 


Calculated  mixing  ratio  of  NaX  (Na  +  NaO  +  NaC>2  +  NaOH  +  NaCl)  for 
different  rates  for  the  polymerization  sink.  The  solid-line  profile 
corresponds  to  the  calculations  in  Fig.  2  (no  polymerization).  The 
dashed  line  assumes  the  polymerization  mechanism  of  Lamb  and  Benson 
(1986)  with  their  estimate  for  the  dimerization  rate  K^.  Cal¬ 
culations  denoted  by  the  dash-dot  line  assume  a  dimerization  rate  a 
factor  of  10  smaller. 


Rates  (cm  3  s  l)  for  conversion  of  HC1  to  active  chlorine  (Cl,  CIO) 
in  the  upper  stratosphere.  Conversion  is  primarily  achieved  in  pre¬ 
sent  models  through  reaction  of  OH  with  HC1  (1).  The  rate  limiting 
step  for  the  sodium-catalyzed  cycle  is  photolysis  of  NaCl  (6).  We 
adopt  the  NaCl  profile  in  Fig.  2  and  the  J  rates  calculated  from 
cross  sections  of  Silver  et  al.  (1986). 


Percent  change  in  calculated  CIO  densities.  The  change  is  calcu¬ 
lated  relative  to  a  baseline  with  no  sodium  chemistry.  Solid, 
dashed,  and  dash-dot  lines  represent  the  same  assumptions  for  the 
sodium  sink  as  in  Fig.  3.  The  open  circles  denote  calculations  with 
twice  the  sodium  abundance  and  no  polymerization  sink. 


Percent  change  in  calculated  O3  densities.  Calculations  correspond 
to  the  CIO  profiles  in  Fig.  3. 


release  of  C?) 


Chemiscry  of  sodium  species  in  the  stratosphere  and  mesosphere. 
Adopted  rates  are  listed  in  Table  1.  Rates  for  processes  denoted 
by  dotted  lines  have  not  been  measured.  Formation  of  (.\'aCl)tl) 
NaHC03 ,  and  attachment  to  aerosols  act  as  permanent  sinks  of 
gaseous  sodium  species. 
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Figure  2.  Calculated  diurnal ly-averaged  densities  of  sodium  species  in  the 
upper  stratosphere.  The  species  Na02  arid  NaCl  disappear  at  night 
in  this  altitude  range,  while  NaCl  exhibits  little  diurnal  varia¬ 
tion.  Loss  of  NaX  is  by  washout  in  the  troposphere.  Other  possi¬ 
ble  sinks  of  NaX  are  neglected.  Model  input  and  assumptions  are 
discussed  in  the  text. 


AO 


ALTITUDE  (KM) 


55k  _ 


—  kd 

--  kd 
•-  kd 


40 


MIXING  RATIO  (pptv) 


Figure  3.  Calculated  mixing  ratio  of  NaX  (Na  +  NaO  +  NaC>2  +  NaOH  +  NaCl)  for 
different  rates  for  the  polymerization  sink.  The  solid-line  pro¬ 
file  corresponds  to  the  calculations  In  Figure  2  (no  polymeriza¬ 
tion).  The  dashed  line  assumes  the  polymerization  mechanism  of 
Lamb  and  Benson  (1986)  with  their  estimate  for  the  dimerization 
rate  K^.  Calculations  denoted  by  the  dash-dot  line  assume  a 
dimerization  rate  a  factor  of  10  smaller  than  in  Lamb  and  Benson 
(  1986). 
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RATE  (cm'V) 

Kates  (cm  3  s  *)  for  conversion  of  HC1  to  active  chlorine  (Cl,  CIO) 
in  the  upper  s t ratosphere .  Conversion  is  primarily  achieved  in 
present  models  through  reaction  of  OH  with  HC1  (1).  The  rate 
limiting  step  for  the  sodium-catalyzed  cycle  is  photolysis  of  NaCl 
( fi ) .  We  adopt  the  NaCl  profile  in  Figure  2  and  the  J  rates 
calculated  from  cross  sections  of  Silver  et  al.  (1986). 
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Figure  5.  Percent  change  in  calculated  CIO  densities.  The  change  is  calcu¬ 
lated  relative  to  a  baseline  with  no  sodium  chemistry.  Solid, 
dashed,  and  dash-dot  lines  represent  the  same  assumptions  for  the 
sodium  sink  as  in  Figure  4.  The  open  circles  denote  calculations 
with  twice  the  sodium  abundance  and  no  polymerization  sink. 
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PART  III 


The  Neutral/Ion  Chemistry  of  Sodium 
In  the  Statosphere  and  Mesosphere 


ABSTRACT 


I 


I 


This  report  presents  results  from  our  model  calculation  of  the  diurnal 
behavior  of  the  neutral  and  ionic  sodium  species.  It  is  shown  that  the 
concentration  of  Na+  above  90  km  exerts  a  direct  influence  on  the  topside 
scale  height  of  the  atomic  sodium  layer  in  the  mesosphere.  The  abundance  of 
Na+  is  controlled  by  the  reaction  rate  for  the  clustering  reaction  which  is 
the  rate  limiting  step  for  the  conversion  of  Na+  to  Na.  Downward  transport  of 
Na+  due  to  ion  drift  is  also  effective  in  reducing  the  sodium  abundance  around 
90  km.  The  model  can  successfully  simulate  the  observed  topside  scale  height 
of  the  sodium  layer  by  either  using  a  slow  rate  for  the  clustering  reaction  or 
incorporating  transport  of  Na+. 
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Figures 


INTRODUCTION 


1  . 

In  Parc  I,  we  presented  results  from  our  diurnal  simulation  of  the  con¬ 
centrations  of  the  sodium  species  based  on  neutral  chemistry.  We  now  extend 
the  discussion  to  include  ion  chemistry  of  sodium  species. 

Pos it ive ly  "charged  ions  of  atomic  weight  23  have  been  observed  in  the 
E  region  during  several  rocket  flights  (see  Swider,  1984,  and  references 
therein).  These  observations  have  been  interpreted  as  positively-charged  Na 
ion  (N'a+).  Sodium  ions  are  produced  by  photo-ionization  as  well  as  by  charge 
exchange  reactions  with  the  major  positive  ions  C>2+  and  N0+  (Swider,  1969a). 
Dissociative  recombination  of  molecular  ions,  produced  by  clustering  of  Na+ 
with  N2  and  CO2  (Richter  and  Sechrist,  1979)  ultimately  removes  Na+  ions.  The 
clustering  rate  decreases  rapidly  with  altitude  resulting  in  inefficient  re¬ 
moval  of  ions  at  high  altitudes.  Although  there  are  insufficient  data  to  ob¬ 
tain  any  detaiLed  information  on  the  diurnal  and  seasonal  behavior  of  the 
ions,  they  can  nonetheless  affect  the  behaviour  of  the  neutral  sodium.  In 
particular,  the  departure  of  the  observed  scale  height  of  atomic  sodium  from 
the  scale  height  of  the  atmosphere  above  90  km  (Gibson  and  Sanford,  1970; 
Simonich  et  ul.,  1979)  is  thought  to  be  controlled  by  the  amount  of  Na+  ions 
present  (Hanson  and  Donaldson,  1967). 

In  order  to  calculate  the  production  and  removal  rates  of  the  sodium 
ions,  one  needs  to  know  the  distributions  of  the  background  ions.  In  Sections 
2  and  3,  we  will  present  the  model  for  calculating  the  diurnal  behavior  of  the 
bickground  Lons  in  the  D  and  lower  E  region.  Section  4  will  present  the  re¬ 
sults  of  the  calculated  diurnal  behavior  of  the  sodium  ions.  The  effects  of 
the  ions  on  the  observations  of  the  scale  height  of  atomic  sodium  above  90  km 
and  other  impact  on  sodium  will  be  discussed  in  Section  5. 

I .  IONS  IN  THE  D  AND  LOWER  E  REGIONS 

Since  our  primary  interest  is  in  the  effect  of  metal  ions'  chemistry  on 
the  neutral  chemistry  around  'll)  km,  we  wi  11  use  a  version  of  the  ion  chemistry 
appropriate  tor  the  i>  u:ni  lower  E  regions.  In  the  lower  E  region,  production 
of  0,  f,  N?+  and  Of  by  x-ray  ant  Lyman-;-!  dominates.  However,  the  N 2+  and  0+ 
are  converted  rapidly  to  Oy  +  :>y  charge  exchange.  Some  0+  is  converted  to  N0+, 
which  is  also  : armed  by  the  charge  exchange  reaction  of  NO  with  0 2+ •  In  the  D 
region,  production  of  V>+  by  Lyman-t  is  important,  making  this  ion  the  most 


important  ionic  constituent  near  80  km.  In  the  lower  D  region  (60-80  km), 
there  are  evidence  from  mass-spectrometric  measurements  that  hydrated  cluster 
ions  are  common  (Narcisi  and  Bailey,  1965;  Goldberg  and  Aikin,  1971;  Krankow- 
sky  et  al.,  1972).  Swider  and  Narcisi  (1983)  suggested  that  carbon  dioxide 
may  be  an  important  agent  in  hastening  the  formation  of  these  clustered  ions. 

Negative  ions  are  produced  through  a  chain  of  reactions  initiated  by 
electron  attachment  to  O2 : 

e  +  O2  +  M  +  O2  +  M.  (2-1) 

Although  observations  show  conflicting  results  making  the  identifications  of 
the  ions  difficult,  it  is  relatively  certain  that  the  major  ions  are  CO3  and 
NO3  (Swider  and  Narcisi,  1983).  For  the  problem  at  hand,  it  is  not  crucial 
that  the  identity  of  the  negative  ions  is  precisely  known. 

Positive  ions  are  destroyed  mainly  by  dissociative  recombination, 

AB+  +  e  •+■  A  +  B ,  (2-2) 

where  AB+  denotes  a  molecular  ion.  Rates  for  dissociative  recombination  are 
usually  very  fast  (~  10  6  cm3  sec"1),  and  result  in  lifetimes  of  the  order  of 
minutes  for  positive  molecular  ions  near  90  km.  Removal  of  negative  ions  is 
effected  by  photodetachraent  reactions 

B~  +  hv  +  B  +  e.  (2-3) 

As  the  concentrations  of  the  negative  ions  become  much  higher  below  80  km  at 
night,  positive-negative  ion  recombination 

A+  +  B  +  products  (2-4) 

becomes  an  important  loss  mechanism. 

A  simplified  model  of  the  reaction  schemes  for  positive  and  negative  ions 
appropriate  for  the  D  and  lower  E  region  are  given  in  Figures  1  and  2  respec¬ 
tively.  The  photo-ionization  sources  incorporated  in  our  model  include  the 


For  daytime  condition: 

a)t  X-rays :  Absorption  and  ionization  cross  sections  of  0,  O2,  and  N2 
are  take  .  from  Swider  (1969b),  for  wavelengths  between  1  and  100A. 
The  solar  flaxes  for  quiet  conditions  are  also  from  Swider  (1969b). 
b,c)  Lyman-g  ionization  of  Op  and  Lyman-q  ionization  of  NO:  Absorption 

cross  sections,  ionization  cross  sections,  and  fluxes  are  taken  from 
Swider  (1969a,  1983)  for  quiet  conditions. 

d)  Energetic  electrons:  Spectra  of  energetic  electrons  at  300  km  mea¬ 
sured  by  the  AE-C  satellite  have  been  presented  by  Torr  et  al. 
(1976).  These  authors  estimated  an  ion  production  rate  of 

10  cm  3  s  *  at  for  nighttime  90  km  although  the  rate  for  daytime 
could  be  considerably  higher.  We  have  assumed  that  the  ion  produc¬ 
tion  profile  has  the  shape  of  the  theoretical  calculation  of  Berger 
et  al.  (1974;  Figure  6,  for  5  kev  e-folding  energy);  this  profile  is 
normalized  to  a  peak  production  rate  of  10  cm  3  s  1  at  90  km.  The 
above  total  production  rate  was  subdivided  into  90 7*  O2  and  \QZ  N0+ 
production  as  suggested  by  Swider  (1969a,  1983)  for  cosmic  rays. 
Because  of  the  large  uncertainties,  production  of  ions  from  ener¬ 
getic  electrons  is  not  included  in  most  of  our  calculations.  How¬ 
ever,  they  are  included  in  Figure  3  for  completeness. 

e)  Ionization  of  O2  (*A):  We  have  used  the  parameterization  suggested 
by  Paulsen  et  al.  (1972). 

f)  Cosmic  rays:  The  parameterization  given  by  Swider  (1983)  for  raid- 
latitude  condition  is  used. 

For  night  time  conditions: 

b,c)  Lyman-8  and  Lyman-q  scattered  by  the  geocorona:  The  scattered 

fluxes  as  a  function  of  altitude  and  solar  zenith  angle  are  taken 
from  Strobel  et  al.  (1974,  1980). 

d)  Energetic  electrons:  see  discussion  for  dayside  ion  condition. 

f)  Cosmic  rays:  As  on  the  dayside. 

A  constant  production  rate  of  10  3  cm  3  s  1  at  all  altitudes  due  to  galactic 
x-rays  (Torkar  and  Friedrich,  1983)  is  included.  The  production  of  ions  by 
starlight  (Strobel  et  al.,  1974,  1980)  has  been  ignored  since  it  is  only  im¬ 
portant  above  ~  110  km. 


^Alphabet  refers  to  labels  on  Figure  3. 
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The  adopted  production  rates  for  ions  are  presented  for  noon  in  Figures 
3(a)  and  3(b)  and  3(c)  for  midnight.  These  calculations  assume  43°  latitude, 
equinox  conditions.  Ionization  by  x-rays  and  Lyman-B  are  dominant  above  an 
altitude  ot  90  km  on  the  dayside.  The  rtightside  sources  are  dominated  by  ion¬ 
ization  by  energetic  electrons  if  values  from  Torr  et  a L .  (1976)  are  adopted. 
We  must  note,  however,  that  this  source  is  quite  variable  (Torr  et  al.,  197b). 
The  profiles  for  ion  production  by  energetic  electrons  presented  in  Figure  3 
should  thus  be  considered  as  upper  bounds.  Values  for  these  production  rates 
smaller  than  ionization  by  Lyman-a  and  Lyman-g  cannot  be  ruled  out.  Uncer¬ 
tainties  also  exist  in  the  production  by  Lyman-a  and  Lyman-B  at  night,  since 
the  fluxes  at  these  wavelengths  are  obtained  from  extrapolation  of  radiative 
transfer  calculations  (Strobel  et  al.,  1974,  1980). 

The  list  of  two-  and  three-body  reactions  included  in  the  model  is  given 
in  Table  1.  The  values  for  the  reaction  rates  are  taken  from  Swider  and  Nar- 
cisi  (1983).  Note  that  we  have  assumed  that  the  production  of  02+  ions  is  due 
only  to  photo-ionization  processes.  The  additional  production  from  charge 
exchange  with  N2  +  and  0+  is  incorporated  into  the  same  term  (Swider,  19b9a). 
The  production  of  N0+  is  from  both  photo-ionizat ion  and  charge  exchange. 


3.  MODEL  FOR  THE  BACKGROUND  IONS 

The  basic  structure  of  our  ion  code  is  Illustrated  by  the  flow  chart  in 
Figure  4.  This  code  was  adapted  from  the  mesopheric  code  described  in  Part 
l.  The  local  concentrat ions  of  the  electron,  the  positive  and  the  negative 
ions,  are  governed  by  the  system  ot  equations: 


[e]  +  [NSUM] 


Li[Ni] 

=  [ PSUM] 


(3-1) 

(3-2) 


where  P| ,  L^ ,  and  [N^]  are  the  local  production  rate,  loss  frequency,  and  num¬ 
ber  density  of  the  i*"*1  ionic  species,  [e]  is  the  electron  density,  and  (NSUM) 
and  [ PSUM J  are  the  sum  of  the  negative  and  positive  ions  concentrat ions ,  re¬ 
spectively.  Equation  (3-1)  corresponds  to  the  local  approximation  of  the  con¬ 
tinuity  equation  for  each  ionic  species  whereas  (3-2)  serves  to  determine  |e) 
by  charge  conservation.  The  numerical  scheme  adopted  to  solve  the  system  can 
he  used  tor  any  arbitrary  grouping  of  the  Ionic  species  and  does  not  assume  a 
chain  structure  tor  the  ion  chemistry. 
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Equation  (3-1)  is  solved  at  a  series  of  discrete  time  steps  using  the 
fully  implicit  scheme  described  by  Sze  et  al.  (1980).  The  initial  values  for 
the  electron  and  positive  ion  densities  at  noon  are  taken  from  our  steady 
state  calculation.  At  other  time  steps,  the  initial  values  are  taken  to  be 
the  electron  and  positive  ion  densities  calculated  at  the  previous  time  step. 
At  each  time  step,  the  equations  are  solved  using  an  iteration  scheme  similar 
to  the  one  used  by  Swider  and  Foley  (1978).  The  values  for  [e]j+i  and 
[PSUMlj+i  for  the  j  +  111^  iteration  are  updated  via  (Swider  &  Foley,  1978) 

[elJ+1  =  {  [  e  ]  j  +'  [PSUM]  j  /(l  +  \)}/2, 

[ PSUM] j+1  =  {(ejj  +  [NSUMjj  +  ( PSUM J j }/2  (3-3) 

with 

X  =  [NSUM] j / [e] j . 

The  process  is  repeated  until  equations  (3-3)  are  satisfied  to  some  specified 
accuracy.  The  solution  is  propagated  forward  using  time  steps  of  one  hour 
during  the  day,  two  hours  at  night,  and  30  minutes  around  sunrise  and  sunset. 
After  24  hours,  the  solutions  are  tested  for  periodicity,  and  the  calculation 
is  repeated  if  periodicity  is  not  achieved  to  within  1%.  Due  to  the  small 
time  constants  for  ionic  processes,  our  solutions  usually  become  periodic  in 
one  to  two  diurnal  cycles. 

The  neutral  background  atmosphere  is  obtained  from  the  diurnal  calcula¬ 
tion  described  in  Part  I.  The  effect  of  the  ion  chemistry  on  neutral 
background  trace  gases  is  negligible  under  most  conditions. 

Results  of  our  calculations  for  the  background  ions  and  electrons  at  45°N 
for  equinox  condition  are  presented  below.  The  diurnal  variation  of  the  elec¬ 
tron  density  is  presented  in  Figure  5  at  altitudes  of  86,  90,  94,  and  100  km. 
Similar  results  for  the  same  altitudes  are  shown  in  Figure  6  for  O2  and  N0+, 
which  are  the  main  positive  ion  constituents  in  this  region.  The  electron 
densities  exhibit  a  day-to-night  ratio  of  about  20  -  30  between  86  and  94  km, 
increasing  to  about  50  at  100  km.  There  is  considerable  variation  during  the 
daytime,  with  the  electron  density  increasing  by  factors  of  2  to  5  between 
sunrise  and  noon.  N0+  densities  are  larger  than  O2  densities  over  most  of 
this  region;  the  diurnal  behavior  of  N0+  is  thus  positively  correlated  with 
that  of  the  electrons.  Below  95  km,  O2  densities  exhibit  a  larger  day-to- 
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night  ratio  (greater  than  200).  This  is  due  to  the  conversion  of  Oj  to  N0+  at 
night  through  charge  exchange  with  NO. 

Our  results  during  the  day  are  in  good  agreement  with  diurnal  calcula¬ 
tions  presented  by  Torkar  and  Friedrich  (1983).  We  note  that  the  magnitude  of 
the  diurnal  variation,  however,  will  depend  on  the  model's  assumptions  for 
nightside  sources.  The  ionization  by  energetic  electrons  exhibits  great  vari¬ 
ability.  In  these  calculations,  we  neglect  ion  production  by  energetic  elec¬ 
trons.  The  Lyman-a  and  Lyman-g  sources  are  computed  from  radiative  transfer 
models  which  become  uncertain  below  ~  100  km  (Strobel  et  al.,  1974,  1980). 

Our  results  should  be  interpreted  only  as  the  general  diurnal  behavior  of 
electrons  and  ions.  More  detailed  information  on  the  nightside  ionization 
sources  at  specific  times  is  needed  to  pinpoint  the  diurnal  variation  to  a 
greater  accuracy. 

Profiles  for  electron  and  total  negative  ion  densities  for  day  and  night 
are  shown  in  Figure  7  for  the  same  conditions  as  in  Figures  5  and  6.  In  light 
of  the  possibly  large  variability  in  the  electron  flux,  we  also  include  the 
calculated  profiles  for  which  the  ionization  from  energetic  electrons  at  night 
is  included  (curves  e,  f).  The  figure  illustrates  that  variability  of  this 
source  can  induce  changes  in  the  electron  density  of  more  than  factors  of  4 
above  90  km. 

The  distributions  of  negative  ions  for  daytime  and  nighttime  conditions 
are  given  in  Figure  8.  The  total  negative  ion  density  in  the  daytime  remains 
smaller  than  the  electron  density  down  to  about  60  km.  During  the  night,  the 
concentration  of  the  negative  ions  becomes  larger  in  the  absence  of  photode¬ 
tachment  reaction  (2-4).  The  density  of  NO3  grows  to  more  than  103  cm  3, 
reflecting  the  increase  in  concentrations  of  NO2  and  O3  which  favors  the 
transfer  to  NO3  via  O3  and  CO3.  Hydrated  negative  ions  are  not  included  in 
the  calculation  since  the  rates  for  the  reactions  involving  these  ions  are 
uncertain. 

Profiles  of  positive  ions  for  day-  and  nighttime  conditions  are  given  in 
Figures  9  and  10.  The  positive  ions  are  mostly  in  the  form  of  O2  and  N0+ 
above  90  km.  Below  90  km,  various  hydrated  ions  are  formed  in  appreciable 
amounts.  Densities  for  these  ions  are  shown  in  Figure  10.  Formation  of  these 
hydrated  ions  is  initiated  by  the  clustering  reactions  (Swider  and  Narclsi, 
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The  profiles  in  Figures  9  and  10  are  calculated  using  the  H2O  mixing  ratios  of 
Bevilacqua  et  al.  (1983),  which  are  up  to  an  order  of  magnitude  smaller  than 
those  of  ALIen  et  al.  (1984)  above  80  km.  Variations  Ln  upper  mesospheric 
water  wiLl  affect  the  calculated  densities  of  the  hydrated  positive  ions.  The 
results  in  Figures  9  and  10  reproduce  the  generaL  characteristics  of  calcula¬ 
tions  by  other  models  (Swider  and  Foley,  1978;  Swider,  1983;  Swider  and  Nar- 
cisi,  1983). 


4.  DIURNAL  MODEL  FOR  SODIUM  NEUTRALS  AND  IONS 

+  +  +  + 

The  ionic  sodium  species  considered  are:  Na  ,  NaO  ,  Na02  ,  Na  '  N2, 

Na+  •  CO2 ,  and  Na+  •  H2O.  The  chemical  scheme  for  this  family  is  illustrated 
in  Figure  11.  This  scheme  combines  the  reactions  proposed  by  Murad  (1978)  and 
Richter  and  Sechrist  (1979).  The  latter  paper  proposed  the  clustering  scheme 
that  recycles  the  Na  +  back  to  neutral  Na. 


The  Na+  Lon  is  produced  by 

Na  +  hv  -*•  Na+  +  e, 

+  + 

Na  +  O2  ■*  Na  +  O2  » 

and 

Na  -t-  No*"  *  Na  +  NO. 

Molecular  ions  are  formed  via 

4"  + 

N  a  +•  Np  +  M  >  Na  •  +  M 

which  in  turn  yields  the  more  stable  Na 
ing  react  ions : 

4-  f- 

Na  •  N2  +•  Ol?  ♦  Na  •  CO?  +  N? 


(4-1) 

(4-2) 

(4-3) 

(4-4) 

•  C(>2  and  Na  •  tl20  through  the  switch- 

(4-5) 
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and 


Na+»  C02  +  H20  +  Na+*  H20  +  C02 .  (4-6) 

The  cluster  ions  are  converted  back,  to  neutral  atomic  sodium  by  dissociative 
recombination.  Reaction  (4-4)  is  the  rate  limiting  step  for  the  sequence  of 
reactions  (4-4)  through  (4-6). 

The  ion  reactions  shown  in  Figure  11  are  added  to  our  diurnal  model  for 
the  sodium  species  (Part  1),  and  the  diurnal  behavior  and  altitude 
distributions  of  the  ionic  and  neutral  sodium  species  are  calculated  in  a 
self-consistent  manner.  In  these  calculations,  the  background  neutral  and 
ionic  constituents  are  not  recalculated,  since  they  are  not  affected  by  the 
sodium  chemistry.  Parameters  used  in  calculating  the  reaction  rates  in  our 
sodium  model  are  shown  in  Table  2,  with  the  rates  for  the  neutral  species  as  ' 
documented  in  Part  I  and  those  for  the  ionic  species  taken  from  Murad  (1978) 
and  Richter  and  Sechrist  (1979).  As  in  the  case  of  the  neutral  chemistry, 
most  of  the  ion  reactions  have  not  been  measured,  and  the  values  proposed  are 
deduced  from  measured  values  of  similar  reactions. 

The  diurnal  variation  of  Na+  is  shown  in  Figure  12  for  90,  94,  and  100 
km.  The  day-to-night  ratio  ranges  from  400  (90  km)  to  8  (100  km).  The  factor 
of  10  variation  for  the  altitude  range  90  to  110  km  suggested  by  the  data  as 
presented  by  Swider  (1984)  is  thus  consistent  with  our  results,  although  the 
large  changes  in  diurnal  variation  over  a  small  altitude  range  preclude  draw¬ 
ing  any  definitive  conclusions  about  the  validity  of  our  chemical  scheme. 

Altitude  profiles  at  noon  for  Na+  and  Na  are  shown  in  Figure  13.  Atomic 
sodium  is  converted  into  Na+  at  high  altitudes  where  their  densities  are  com¬ 
parable.  As  we  shall  see  in  the  next  section,  processes  that  affect  the  con¬ 
centration  of  Na+  can  then  indirectly  alter  the  distributions  of  Na  at  these 
altitudes.  Densities  of  the  other  ionic  species  are  smaller  by  factors  of  100 
or  more.  In  particular,  the  [NaO+]/[Na+]  value  is  ~  10  7  to  10  9,  in  agree¬ 
ment  with  the  estimations  of  Murad  (1978). 

5.  SODIUM  ION  CHEMISTRY  AND  THE  SCALE  HEIGHT  OF  THE 

NEUTRAL  ATOMIC  SODIUM  LAYER 

The  topside  scale  height  of  the  sodium  layer  has  been  observed  to  be  much 
smaller  than  that  of  the  background  atmosphere  (Gibson  and  Sandford,  1971; 
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Megie  and  Blaraonc,  1977;  Simonich  et  al.,  1979;  Granier  and  Megie,  1982),  de¬ 
creasing  Co  abouC  2-4  km  at  100  km  altitude.  A  possible  explanation  for  this 
small  scale  height  is  the  conversion  of  Na  to  Na+  by  photoionizat Lon  or  charge 
exchange  with  02 +  (Hanson  and  Donaldson,  1967;  Gadsden,  1968;  Richter  and 
Sechrist,  1979). 

This  sensitivity  of  the  scale  height  of  atomic  Na  ion  chemistry  is  illus¬ 
trated  in  Figure  14  which  gives  the  calculated  profiles  of  Na  with  different 
assumptions  for  the  Na+  abundance.  As  discussed  in  Section  4,  reaction  (4-4) 
constitutes  the  rate-limiting  step  for  recombination  of  sodium  ions  and  is 
assumed  to  be  of  the  form  (Richter  and  Sechrist,  1979): 


;<s  =  K*  cm*  ,-l 


(5-1) 


A  value  of  K*  =  2  x  IO-31  cm6  s_1  was  suggested  by  Richter  and  Sechrist  (1979) 
by  analogy  to  similar  reactions  involving  N0+.  Reaction  (4-4)  has  also  been 
measured  at  110°  K  in  a  helium  buffer  by  Perry  et  al.  (1980);  their  value  of 
2.6  x  10  30  cm0  s  1  implies  a  value  of  K*  =  7  x  10  31  cm6  s_1  if  the  tempera¬ 
ture  dependence  in  eq .  (5-1)  is  valid.  The  Na  profiles  in  Figure  14  are  cal¬ 
culated  using  different  values  for  K*.  The  noontime  profile  from  the  calcula¬ 
tion  without  ionic  species  (-0-)  is  compared  with  the  cases  calculated  with 
value  of  K*  from  Richter  and  Sechrist  (1979)  (K*  =  2  x  10"31  cm6  s  )  and 

Perry  et  al.  (1980)  (K*  =  7  x  10  3'  cm*3  s  3,  - ).  We  have  also  considered 

the  hypothetical  situation  with  K*  =  4  x  10  32  cm6  s_1  (9G9i) .  As  expected, 
the  scale  height  of  neutral  sodium  increases  with  increasing  K* ,  reflecting 
the  more  efficient  rate  for  reconversion  of  Na+  to  Na.  The  observed  small 
scale  heights  about  90  km  could  thus  be  explained  by  the  ion  chemistry  without 
having  to  invoke  transport  processes  (as  in  Hanson  and  Donaldson,  1967,  for 
example)  if  Kg  is  sufficiently  fast. 

Local  values  for  the  topside  scale  height  are  given  in  Table  3  for  the 
cases  considered  in  Figure  14.  A  comparison  with  observations  indicates  an 
upper  limit  of  2  x  10  31  ciif’  s  *  for  K*  if  the  small  topside  scale  heights  are 
to  be  explained  solely  by  ion-neutral  chemistry.  Such  small  values  of  K*  are 
not  likely,  however,  given  the  measurements  of  Perry  et  al.  (1980)  and  that 
the  faster  clustering  rate  are  expected  when  N2  is  the  buffer  gas  (Jegou  et 
al.,  1984;  Heimerl,  197t>).  An  alternative  mechanism  to  reduce  the  topside 
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scale  height  Is  the  downward  transport  of  Na+  from  the  region  prior  to  its 
conversion  to  Na.  Vertical  ion  drifts  with  velocities  of  10  to  50  era  s_1  have 
been  suggested  by  Hunten  (1981)  and  Kirchoff  and  Batista  (1984)  and  could  be 
effective. 

Ionic  species  could  acquire  a  vertical  velocity  through  electromagnetic 
interactions.  These  include  electrostatic  forces  from  the  ionospheric  elec¬ 
tric  fields  (Kirchoff  and  Batista,  1984)  and  v  x  B  type  forces  from  the  mag¬ 
netic  field  because  of  the  motions  in  the  horizontal  plane  associated  with  the 
zonal  wind  (Jegou  et  al.,  1984)  and  tides  (Kirchoff  and  Batista,  1984).  Ion 
drift  could  provide  an  effective  way  for  the  vertical  transport  of  Na  species 
in  the  mesosphere,  particularly  when  a  significant  portion  of  the  sodium  spe¬ 
cies  is  in  ionic  form. 

A  detailed  calculation  of  these  effects  would  require  the  incorporation 
of  the  above  parameters  in  our  model.  In  lieu  of  a  full  calculation,  we  exam¬ 
ined  the  response  of  the  topside  layer  to  typical  net  vertical  downward  velo¬ 
cities  for  ions  by  modifying  the  continuity  equation  for  the  NaX  species  to 
have  the  form 


i!  +  d[  *  iNax]  •  wd]  . 


P  -  L, 


(5-2) 


where  $  is  the  flux  due  to  molecular  and  eddy  diffusion,  P  and  L  denote  the 
production  and  loss  mechanisms  for  NaX,  and  represents  the  vertical  ion 
drift  velocity.  The  relative  abundance  [Na+]/[NaX]  is  independent  of  the  as¬ 
sumed  drift  velocity  if  the  time  constant  for  the  vertical  drift  is  larger 
than  the  time  constant  for  the  removal  of  Na+  as  given  by  the  rate  limiting 
step  (4-4).  This  condition  is  valid  below  100  km,  where  the  time  constant  for 
the  vertical  drift  is  ~  101*  s,  while  the  time  constant  for  Na+  clustering  is 
~  103  s. 

Figure  15  shows  the  results  of  our  calculations  for  constant  drift  velo¬ 
cities  equal  to  0,  -5,  -10,  and  -20  cm/s  between  120  and  90  km.  The  small 
topside  scale  height  observed  (2-3  km,  Gibson  and  Sandford,  1981;  Cleraesha  et 
al.,  1982)  can  be  modeled  using  downward  velocities  between  10  and  20  cm/s. 

The  effect  of  vertical  ton  drift  is  thus  an  important  contributing  factor  in 
controlling  the  topside  sodium  layer. 
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6 .  CONCLUDING  REMARKS 

We  have  presented  the  results  of  a  fully  diurnal  model  for  ion  species  in 
the  D  and  lower  E  regions.  When  coupled  with  our  previous  model  of  the  neu¬ 
tral  species,  our  present  model  allows  for  the  self-consistent  calculation  of 
background  neutral  and  ionic  species  between  60  and  120  km.  Results  of  our 
calculations  are  in  good  agreement  with  those  of  other  models. 

Production  of  Na+  occurs  primarily  through  photoionization  below  iOO  km, 
and  through  charge  exchange  with  02+  at  higher  altitudes.  The  main  loss 
mechanism  is  dissociative  recombination  of  sodium  cluster  ions.  A  self- 
consistent  calculation  thus  requires  knowledge  of  the  02+  and  e  densities  in 
the  above  altitude  regimes.  Our  code  has  been  applied  to  the  calculation  of 
the  density  and  diurnal  variation  of  sodium  species  above  90  km.  The  calcu¬ 
lated  diurnal  variation  of  Na+  (factor  of  10-100  between  90  and  100  km)  and 
daytime  densities  (~  102  cm  3 )  are  in  good  agreement  with  rocket  data  (Swider, 
1984). 

We  have  also  investigated  the  effect  of  ion  chemistry  on  the  abundance  of 
neutral  sodium  above  90  km.  In  particular,  the  small  scale  height  (2-3  km) 
observed  in  the  topside  sodium  layer  seems  to  require  either  downward  ion 
drifts  velocities  of  order  10-20  cm  s  1  above  90  km  or  a  rate  of  clustering  of 
Na+  with  N2  a  factor  of  5  smaller  than  the  value  suggested  by  Richter  and 
Sechrist  (1979).  Further  elucidation  of  this  question  will  require  measure¬ 
ments  of  the  clustering  rate  at  mesospheric  temperatures  as  well  as  incorpora¬ 
tion  of  more  realistic  scenarios  for  ion  drifts  into  our  model. 

Since  we  have  included  a  complete  treatment  of  the  ion  chemistry  in  the  D 
region,  the  model  presented  above  is  also  suited  to  studies  of  the  D  region 
ion  chemistry  and  its  possible  effect  on  neutral  constituents  during  disturbed 
conditions  (Swider  and  Keneshea,  1973;  Solomon  et  al.,  1983).  The  model  also 
allows  for  each  extension  to  the  calculation  of  other  metallic  ion  species. 
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List  of  Reactions  and  Reaction  Rates  tor  the  Background  Ions 


TYPE3 

PARI 

PAR2 

PAR3 

REACT  I  UN'S 

1 

7. 

l.OOE-28 

-3. OOE+OO 

O.OOE+OO 

H+.(H20)3  +  CU2  +  M  =  lif.  (  I2U)  3.C02  +  M, 

2 

1. 

1.00E-09 

O.OOE+OO 

O.OOE+OO 

H+.(H20)3.C02  +  H20  =  I1+.(H20)4  +  C02,  1 

3 

1. 

4.50E-10 

0.00E+00 

O.OOE+OO 

02+  +  NO  =  N0+  +  02,  1 

4 

7. 

1.60E-30 

-3.20E+00 

O.OOE+OO 

02  +  02  +  M  =  O4 +  +  M ,  7 

5 

1. 

3.00E-10 

0.00E+00 

O.OOE+OO 

04  +  0  =  02+  +  03 ,  1 

6 

1 . 

1.30E-09 

0.00E+00 

O.OOE+OO 

o4+  +  h2o  =  02+.II2()  +  02,  1 

7 

1. 

1.00E-09 

O.OOE+OO 

O.OOE+OO 

02+.H20  +  H20  =  H+.H200H  +  02,  1 

8 

1. 

2 . 00E-10 

0.00E+00 

O.OOE+OO 

02+.H20  +  H20  =  H+.H20  +  OH,  1 

9 

1. 

1.40E-09 

0.00E+00 

O.OOE+OO 

H+.H2OOH  +  H20  =  H+.(H20)2  +  OH,  1 

0 

7. 

L.40E-30 

-2.00E+00 

O.OOE+OO 

H+.H20  +  n2  +  n2  =  H+.H2O.N2  +  N2,  7 

l 

1. 

1.00E-09 

O.OOE+OO 

O.OOE+OO 

H+.n2O.N2  +  co2  =  u+.h2o.co2  +  n2,  1 

2 

1. 

1.00E-09 

0.00E+00 

O.OOE+OO 

H+.H20.C02  +  H20  =  H+.(H20)2  +  co2,  1 

3 

7. 

5.00E-28 

-3.00E+00 

O.OOE+OO 

H+.(H20)2  +  C02  +  M  =  H+.(H20)2C02  +  M,  7 

4 

1. 

1.00E-09 

O.OOE+OO 

O.OOE+OO 

H+.(H20)2C02  +  H20  =  H+.(H20)3  +  co2,  1 

5 

7. 

3.50E-27 

-3.00E+00 

O.OOE+OO 

H+.H20  +  H20  +  M  =  H+.(H20)2  +  M,  7 

b 

7. 

2.20E-27 

-3.00E+00 

0. OOE+OO 

H  .(H20)2  +  H20  +  M  =  H+.(H20)3  +  M,  7 

7 

7. 

2.30E-27 

-3.00E+00 

O.OOE+OO 

H+.(H20)3  +  H20  +  M  =  H+.(H20)4  +  M,  7 

8 

7. 

1.00E-29 

-2. 30E+00 

O.OOE+OO 

N0+  +  C02  +  N2  =  N0+.C02  +  N2,  7 

9 

1. 

1.00E-09 

O.OOE+OO 

O.OOE+OO 

no+.co2  +  H20  =  C02  +  N0+.H20,  1 

:u 

7. 

l . 10E-27 

-4. 70K+00 

O.OOE+OO 

no+.h2o  +  h2o  +  n2  =  N0+.(H20)2  +  N2,  7 

:i 

7. 

L.60E-27 

-4.70E+00 

O.OOE+OO 

no+.(h2o)2  +  h2o  +  n2  =  no+.(h2o)3  +  n2. 

2 

1. 

7.00E-11 

O.OOE+OO 

O.OOE+OO 

no+.(h2o)3  +  h2o  =  H+.(H20)3  +  XX,  I 

3 

7. 

4.00E-07 

-1  .OOE+OO 

O.OOE+OO 

N0+  +  E  =  N  +  0,  7 

:4 

7. 

2. L0E-07 

-6.30E-01 

O.OOE+OO 

02+  +  E  =  0  +  0,  7 

:5 

1. 

2.00E-06 

O.OOE+OO 

O.OOE+OO 

04+  +  E  =  02  +  02,  l 

aThe  parameterization  and  values  for  the  reaction  constants  are  taken  from 
Swlder  and  Narcisi  (1983).  The  reaction  rates  are  calculated  as  follows: 

Type  l:  k  =  FAR  1 

Type  7:  k  =  PAK1*( 300/T)_PAK2 

Type  8:  k  =  PAR  I  *(  300/T )-l>AK2*exp(  PAR 3/T  ) . 


type 

PARI 

PAR2 

PAR3 

REACTIONS 

26 

1. 

1.50E-06 

0.00E+00 

0. OOE+OO 

02+-H20  +  E  =  O2  +  H2O,  1 

27 

1 . 

2.00E-06 

0.00E+00 

0. OOE+OO 

H+.H200H  +  E  =  H2O  +  H2O,  1 

28 

1. 

1.30E-06 

0. 00E+00 

0. OOE+OO 

H+.H20  +  E  =  H2  +  OH,  1 

29 

1. 

2.50E-06 

0 . OOE+OO 

0. OOE+OO 

H+.(H20)2  +  E  =  H20  +  H2  +  OH,  1 

30 

1 . 

3.00E-06 

0.00E+00 

0. OOE+OO 

H+ .  ( H20 )  3  +  E  =  OH  +  XX,  1 

31 

1 . 

3.60E-06 

0. OOE+OO 

0. OOE+OO 

H+.(H20)4  +  E  =  OH  +  XX,  1 

32 

1. 

1 .50E-06 

0.  OOE+OO 

0. OOE+OO 

H+.H2O.N2  +  E  =  H  +  XX,  1 

33 

1 . 

2.00E-06 

0. OOE+OO 

0. OOE+OO 

H+.H20.C02  +  E  =  H  +  XX,  1 

34 

1 . 

3.00E-06 

0. OOE+OO 

0. OOE+OO 

H+.(H20)2C02  +  E  =  H  +  XX,  1 

35 

1 . 

1.50E-06 

0. OOE+OO 

0. OOE+OO 

N0+.C02  +  E  =  NO  +  co2,  1 

36 

1 . 

l .50E-06 

0. OOE+OO 

0. OOE+OO 

NO+.H2C  +  E  =  NO  +  H20,  1 

37 

1. 

2.00E-06 

0. OOE+OO 

0. OOE+OO 

N0+.(H20)2  +  E  =  NO  +  XX,  1 

38 

1. 

3.00E-06 

0. OOE+OO 

0.  OOE+OO 

N0+.(H20)3  +  E  =  NO  +  XX,  1 

39 

8. 

3.30E-06 

-4.  OOE+OO 

-5.03E+03 

04+  +  M  =  02+  +  02  +  M,  8 

40 

7. 

2.00E-31 

-4.40E+00 

0. OOE+OO 

N0+  +  N2  +  N2  =  N0+.N2  +  N2,  7 

41 

1. 

l .00E-09 

0. OOE+OO 

0. OOE+OO 

no+.n2  +  co2  =  no+.co2  +  n2,  1 

42 

8. 

1 .00E-08 

-5.40E+00 

-2. 10E+03 

N0+.N2  +  N2  =  N0+  +  N2  +  N2,  8 

43 

7. 

2.00E-31 

-4.40E+00 

0.  OOE+OO 

no+.h2o  +  n2  +  n2  =  NO+.H2O.N2  +  N2,  7 

44 

1. 

1 .00E-O9 

0. OOE+OO 

0. OOE+OO 

NO+.H2O.N2  +  C02  =  N0+.H20C02  +  N2,  1 

45 

1 . 

1 .00E-09 

0. OOE+OO 

0. OOE+OO 

N0+.H20.C02  +  H20  =  N0+.(H20)2  +  co2,  1 

46 

7. 

2.00E-31 

-4.40E+00 

0. OOE+OO 

N0+ . ( H20 ) 2  +  N2  +  N  2=  N0+.(H20) 2.N2  +  N2 

47 

1 . 

1.00E-09 

0. OOE+OO 

0. OOE+OO 

N0+.(H20)2.N2  +  C02=  N0+.(H20) 2.C02  +  n2 

48 

1 . 

l .00E-09 

0. OOE+OO 

0.  OOE+OO 

no+.(h2o)2.co2  +  H20  =  N0+.(H20)3  +  C02, 

49 

8. 

1.20E-0I 

-4. OOE+OO 

-8.80E+03 

H+.(H20)4  +  M  =  H+.(H20)3  +  H20  +  M,  8 

50 

8. 

2.00E-05 

-3.30E+00 

-5. 50E+03 

N0+.C02  +  M  =  N0+  +  C02  +  M,  8 

51 

8. 

5.00E-08 

-5.40E+00 

-2. 10E+03 

N0+.H20.N2  +  M  =  N0+.H20  +  n2  +  M,  8 

52 

8. 

2.00E-04 

-3.30E+00 

-5. 50E+03 

N0+.H20.C02  +  M  =  NO+.H20  +  C02  +  M,  8 

3  3 

8. 

5.0OE-07 

-5.40E+00 

-2. 10E+03 

N0+.(H20)2.N2  +  M  =  N0+.(H20)2  +  N 2  +  M, 

34 

8. 

2.00E-03 

-3.30E+00 

-5.50E+03 

N0+ . ( H20 ) 2 • C02+  M  =  N0+.(H20)2  +  C02  +  M 

5  3 

8. 

1 . OOE-02 

-5.70E+00 

-6. 8 OK +0  3 

NU+.(H20)3  +  M  =  N0+.(H20)2  +  H20  +  M,  8 

56 

8. 

1 .00E-02 

-5.70E+00 

-8. 10E+03 

N0+ . ( H20 ) 2  +  M  =  N0+.H20  +  H^  +  M,  8 

57 

1 . 

5.00E-10 

0.  OOE+OO 

0.  OOE+OO 

02  +03=03  +02,  1 

38 

8. 

1 . 40E-29 

-1  .OOE+OO 

-8.00E+02 

02  +  E  +  0 2  =  02  +  02,  8 

59 

1 . 

1.00E-31 

0.  OOE+OO 

0.  OOE+OO 

02+E+N2=02  +  N  2  ,  1 

58 


59  1  . 


TYPE  PARI  PAR2  PAR3  REACTION 

60  1.  I.50E-10  O.OOE+OO  0.00E+00  02“  +  0  =  O3  +  E»  1 

61  1.  1.50E-10  O.OOE+OO  O.OOE+OO  02"  +  0  =  0~  +  02 ,  1 

62  l.  2.00E-10  O.OOE+OO  O.OOE+OO  02~  +  02(1A)  =  02  +  02  +  E,  1 

63  7.  2.50E-30  -2.50E+00  O.OOE+OO  02“  +  02  +  02  =  04~  +  02,  7 

63  l.  4.30E-10  O.OOE+OO  O.OOE+OO  04“  +  C02  =  C04~  +  02,  1 

65  1.  1.30E-I0  O.OOE+OO  O.OOE+OO  C04~  +  03  =  03“  f  C02  +  02)  I 

66  1.  4.80E-11  O.OOE+OO  O.OOE+OO  C04~  +  NO  =  N03~  +  C02,  1 

67  I.  1.50E-10  O.OOE+OO  O.OOE+OO  C04-  +  0  =  C03"  +  02,  1 

68  I.  3.20E-10  O.OOE+OO  O.OOE+OO  03_  +  0  =  02“  +  02,  1 

69  1.  5.50E-10  O.OOE+OO  O.OOE+OO  03“  +  C02  =  C03~  +  02,  1 

70  1.  l.LOE-LO  O.OOE+OO  O.OOE+OO  CO3-  +  0  =  02“  +  C02,  1 

71  1.  1.10E-11  O.OOE+OO  O.OOE+OO  Co3-  +  NO  =  N02~  +  C02,  1 

72  L.  2.00E-L0  O.OOE+OO  O.OOE+OO  0)3"  +  N02  =  NO3-  +  C02,  l 

73  1.  9.00E-11  O.OOE+OO  O.OOE+OO  N02“  +  03  =  N03~  +  02>  1 

74  7.  5.00E-30  -1.90E+00  O.OOE+OO  O-  +  02  +  02  =  03"  +  02,  7 

75  1.  2.00E-I0  O.OOE+OO  O.OOE+OO  0~  +  0  =  02  +  E,  l 

76  l.  3.00E-I0  O.OOE+OO  O.OOE+OO  0~  +  02(lA)  =  03  +  E,  I 

77  1.  4.40E-10  O.OOE+OO  O.OOE+OO  0“  +  03  =  03~  +  0,  1 

78  1.  4.40E-10  O.OOE+OO  O.OOE+OO  0~  +  03  =  02“  +  02,  1 

79  1.  6.00E-08  O.OOE+OO  O.OOE+OO  02+  +  [NSUM]b  =  XX,  1 

80  l.  6.00E-08  O.OOE+OO  O.OOE+OO  04+  +  [NSUM]  =  XX,  1 

81  I.  6.00E-08  O.OOE+OO  O.OOE+OO  02  +  H20  +  [NSUM]  =  XX,  l 

82  1.  6.00E-08  O.OOE+OO  O.OOE+OO  H+.H2O.OH  +  [NSUM]  =  XX,  1 

83  l.  6.00E-08  O.OOE+OO  O.OOE+OO  H+.H20  +  [NSUM]  =  XX,  l 

84  1.  6.00E-08  O.OOE+OO  O.OOE+OO  H+.H2O.N2  +  [NSUM]  =  XX,  l 

85  1.  6.00E-08  O.OOE+OO  O.OOE+OO  H+.H20.C02  +  [NSUM]  =  XX,  l 

86  1.  6.00E-08  O.OOE+OO  O.OOE+OO  H.(H2Q)2+  +  [NSUM]  =  XX,  I 

87  1.  6.00E-08  O.OOE+OO  O.OOE+OO  H.(H20)2.C02  +  [NSUM]  =  XX,  l 

88  1.  6.00E-08  O.OOE+OO  O.OOE+OO  H+.(H20)3  +  [NSUM]  =  XX,  1 

89  1.  6.00E-08  O.OOE+OO  O.OOE+OO  H+. ( H20) 3 .C02  +  [NSUM]  =  XX,  1 

90  1.  6.00E-08  O.OOE+OO  O.OOE+OO  H+.(H20)4  +  [NSUM]  =  XX,  L 

91  1.  6.00E-08  O.OOE+OO  O.OOE+OO  N0+  +  [NSUM]  =  XX,  l 

b[NSl'M]  denotes  the  total  density  of  negative  Ions. 
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IV  Kb 

PARI 

PAR2 

PAR  f 

R8ACTL0N 

1  . 

6.008-08 

0.008+00 

0.008+00 

N0+.N2  +  [NSUM]  =  XX, 

1 

i  . 

6.008-08 

0.008+00 

0.008+00 

N0+.C02  +  [NSUM]  =  XX, 

1 

1. 

6 .008-08 

0.008+00 

0.008+00 

02~  +  [ PSUM] c  =  XX,  1 

1  . 

6.008-08 

0.008+00 

0.008+00 

0“  +  [PSUM]  =  XX,  1 

1  . 

6.008-08 

0.008+00 

0.008+00 

04"  +  [PSUM]  =  XX,  1 

i. 

6.008-08 

0.008+00 

0.008+00 

03~  +  [PSUM]  =  XX,  l 

1  . 

6.008-08 

0.008+00 

0.008+00 

C03~  +  [PSUM]  =  XX,  l 

1  . 

6.008-08 

0.008+00 

0.008+00 

C04~  +  [PSUM]  =  XX,  1 

L . 

6.008-08 

0.008+00 

0.008+00 

N02~  +  [ PSUM]  =  XX,  1 

1 . 

6.008-08 

0.008+00 

0.008+00 

NO 3"  +  [PSUM]  =  XX,  1 

1 . 

6.008-08 

0.008+00 

0.008+00 

N0+.H20  +  [NSUM]  =  XX, 

1 

1 . 

6.008-08 

0.008+00 

0.008+00 

N0+.(H20)2  +  [NSUM]  = 

XX, 

L 

1  . 

6.008-08 

0.008+00 

0.008+00 

N0+.(H20)3  +  [NSUM]  = 

XX, 

1 

1 . 

6.008-08 

0.008+00 

0.008+00 

NO+.H2O.N2  +  [NSUM]  = 

XX, 

1 

1 . 

6.008-08 

0.008+00 

0.008+00 

N0+.H20.C02  +  [NSUM]  = 

XX 

,  1 

1 . 

6.008-08 

0.008+00 

0.008+00 

N0+.(H20)2 .N2  +  [NSUM] 

= 

XX,  1 

1 . 

6.008-08 

0.008+00 

0.008+00 

N0+.(H20)2.C02  +  [NSUM]  = 

XX,  1 

[PSUM]  denotes  the  total  ilensity  ol  positive  ions. 
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Table  2 

List  of  React  loos  for  the  Sodium  Species 


type3 

PARI 

PAR2 

PAR3 

REACTIONS 

1 

1. 

3.00E-10 

0.00E+00 

0.00E+00 

Na  +  03  = 

NaO  +02,  1 

2 

1. 

4.00E-11 

0.00E+00 

0.00E+00 

NaO  +  0  = 

Na  +  02>  1 

3 

2. 

8.00E-31 

-2.90E+02 

0.00E+00 

Na+02  +  M 

=  NaO 2  +  M,  2 

4 

1. 

1.00E-13 

0.00E+00 

0.00E+00 

Na02  +  0  = 

NaO  +  02»  1 

5 

1. 

0.0QE+00 

0.00E+00 

0.00E+00 

NaO  +  H20 

=  NaOH  +  OH,  1 

6 

1. 

3.00E-11 

0.00E+00 

O.OOE+OO 

NaO  +  03  = 

Na02  +  02»  1 

7 

1. 

1.00E-14 

0.00E+00 

0.00E+00 

NaO  +  H  = 

Na  +  OH,  1 

8 

1. 

1 .00E-11 

0.00E+00 

O.OOE+OO 

Na02  +  OH 

=  NaOH  +  02 >  1 

9 

1. 

1.00E-13 

0.00E+00 

0.00E+00 

Na02  +  H  = 

NaOH  +  0,  1 

10 

1. 

1.40E-12 

0.00E+00 

0.00E+00 

NaOH  +  H  = 

Na  +  H2O,  1 

11 

1. 

1.00E-10 

0.00E+00 

O.OOE+OO 

NaOH  +  0*D 

=  NaO  +  OH,  1 

12 

1. 

1.40E-11 

0.00E+00 

O.OOE+OO 

Na  +  H02  = 

NaOH  +  0,  1 

13 

1. 

1.00E-10 

0.00E+00 

0.00E+00 

NaO  +  03  = 

Na  +  XX,  l 

14 

1. 

1.00E-11 

0.00E+00 

0.00E+00 

NaO  +  H02 

=  NaOH  +  XX,  1 

13 

1. 

6.00E-10 

0.00E+00 

0.00E+00 

Na  +  02+  = 

Na+  +  02,  1 

16 

1. 

8.00E-1 1 

0.00E+00 

0.00E+00 

Na  +  N0+  = 

Na+  +  NO,  1 

17 

7. 

2.00E-31 

-1.50E+00 

0.00E+00 

Na+  +  N2  + 

M  =  Na+.N2  + 

M,  7 

18 

8. 

9.00E-07 

-2.50E+00 

-4.75E+03 

Na+.N2  +  M 

=  Na+  +  N2  + 

N  2>  8 

19 

7. 

2.00E-29 

-1 .50E+00 

O.OOE+OO 

Na+  +  C02 

+  M  =  Na+.C02 

+  M,  7 

20 

8. 

6.00E-06 

-2.50E+00 

-6.25E+03 

Na+ . C02  + 

M  =  Na+  +  C02 

+  M,  8 

21 

7. 

5.00E-32 

-1.50E+00 

O.OOE+OO 

Na+  +  02  + 

M  =  Na02+  +  M 

,  7 

22 

8. 

8.00E-09 

-2.50E+00 

-2.85E+03 

Na02+  +  M 

=  Na+  +  02  +  M 

,  « 

23 

7. 

1.00E-28 

-1.50E+00 

O.OOE+OO 

Na+  +  1120 

+  M  =  Na+.H20 

+  M,  7 

aThe  parameterization  and  values  for  the  reaction  constants  are  taken  from 
Murad  (1978)  and  Richter  and  Schrist  (1979).  The  reaction  rates  are 

calculated  as  follows: 


Type 

1 : 

k  =  PARI 

Type 

2: 

k  =  PARl*exp(-PAR2/T ) 

Type 

7  : 

k  =  PAR1*(300/T)-PAR2 

Type 

8: 

k  =  PARl*(300/T)-PAR2*exp(PAR3/T) 

TYPE 

PARI 

PAR  2 

PAR3 

REACTIONS 

1. 

1.00E-09 

O.OOE+OO 

O.OOE+OO 

Na+.N2  +  C02  =  Na+. 

CO2  +  N 

1. 

I.00E-09 

0.00E+00 

O.OOE+OO 

Na+.C02  +  H20  =  Na+ 

.H20  +  l 

7. 

L.00E-06 

-5.00E-01 

O.OOE+OO 

Na+.C02  +  E  =  Na  + 

co2,  7 

7. 

1.00E-06 

-5.00E-01 

O.OOE+OO 

Na+.H20  +  E  =  Na  + 

H20,  7 

7. 

1 .00E-06 

-5.00E-01 

O.OOE+OO 

Na^.Nj  +  E  =  Na  +  N 

2,  7 

7. 

I.00E-06 

-5.00E-01 

O.OOE+OO 

Na02+  +  E  =  Na02  + 

,  7 

1. 

1 .OOE-IO 

O.OOE+OO 

O.OOE+OO 

Na02+  +  0  =  NaO+  + 

o2,  1 

1. 

5.00E-10 

O.OOE+OO 

O.OOE+OO 

NaO+  +  0  =  Na+  +  02 

.  1 

1. 

1 .00E-U9 

O.OOE+OO 

O.OOE+OO 

NaO  +  02+  =  NaO4  + 

o2,  1 

1. 

1.00E-09 

O.OOE+OO 

O.OOE+OO 

NaO  +  NO+  =  NaO+  + 

NO,  1 

1. 

4  .(JOE-07 

O.OOE+OO 

O.OOE+OO 

NaO+  +  E  =  Na  +  0, 

1 

1. 

6.00E-08 

0.00E+00 

O.OOE+OO 

Na+  +  [ NSUM] b  =  Na 

+  XX,  1 

1. 

6.00E-08 

O.OOE+OO 

O.OOE+OO 

Na+ . N2  +  [NSUM]  =  Na  +  XX, 

1. 

6.00E-U8 

0.00E+00 

O.OOE+OO 

Na+ ,C02  +  [NSUM]  = 

Na  +  XX 

1. 

6.00E-08 

O.OOE+OO 

O.OOE+OO 

Na+.H20  +  [NSUM]  = 

Na  +  XX 

1. 

6.00E-08 

O.OOE+OO 

O.OOE+OO 

NaO+  +  [NSUM]  =  NaO 

+  XX, 

1. 

6.00E-08 

O.OOE+OO 

O.OOE+OO 

Na02+  +  [NSUM]  =  Na 

+  XX, 

NSUM|  denotes  the  total  density  of  negative  ions. 


Table  3 

Calculated  Local  Topside  Scale  Heights  (km) 


1 

Altitude  | 

Scale  Height^ 

(km) 

(km)  | 

No  ions 

K*  =  7  x  10-3ltt  K*  = 

2  x  10~  3 1 

K*  =  4  x  10  32 

1 

92  | 

9.2 

8.13 

7.84 

6.85 

94  | 

6.72 

6.52 

6.11 

5.23 

96  | 

6.13 

5.74 

5.12 

4.31 

98  | 

5.85 

4.99 

4.27 

3.31 

100  | 

_ L 

5.73 

4.26 

3.60 

2.48 

TLocal  scale  height  is  obtained  from  the  gradient  of  Na  profile. 

tTNa+  +  N2  +  M  +  Na+  ♦  N2  +  M  with  reaction  rate  Kg  =  K*  •  Units  for 

K*  are  cm6  s  1  . 

Experimental  Results 

Gibson  and  Sandford  (1970):  4  km  (98  km),  3  km  (100  km),  2  km  (102  km). 

Megie  and  Blamont  (1977):  2-5  km  average  scale  height  above  peak. 

Simonich  et  al.  (1979):  2.5-5  km  at  100  km  (scale  height  changes  during 
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Figure  Captions 


Figure  1.  Schematic  diagram  for  the  chemical  scheme  of  the  positive  ions  in 
the  D  and  lower  E  regions.  The  ions  are  grouped  into  two  families 
with  O2 +  and  N0+  being  the  primary  ions  in  each  group. 

Figure  2.  Schematic  diagram  for  the  chemical  scheme  of  the  negative  ions  in 
the  D  and  lower  E  regions. 

Figure  3.  The  production  rates  for  positive  ions  for  daytime  and  nighttime 

condlt Lons :  A)  Day  (45°  Solar  Zenith  Angle  (SZA)),  altitude  range 
90-120  km);  B)  Day  (45°  SZA),  altitude  range  60-90  km;  c)  Night 
(180°  SZA),  altitude  range  60-120  km.  An  altitude-independent 
production  rate  of  10  3  cra3  s  1  ^ue  to  gallactic  x-rays  is  also 
included  in  the  model.  The  curves  are:  a)  x-ray,  b)  Lyman-B,  c) 
Lyman-a  ,  d)  energetic  electrons,  e)  ionization  of  O2OA),  and  f) 
cosmic  rays. 

Figure  4.  Flow  chart  for  the  diurnal  model  of  neutral  and  ionic  sodium 
species. 

Figure  5.  Calculated  diurnal  variation  of  electron  density.  Profiles  are 
calculated  at  (a)  86  km,  (b)  90  km,  (c)  94  km,  and  (d)  100  km. 

Figure  6.  Calculated  diurnal  variation  of  the  densities  of  NO+  and  02^. 

Profiles  are  calculated  at  (a)  86  km,  (b)  90  km,  (c)  94  km,  and 
(d)  100  km. 

Figure  7.  Calculated  total  electron  and  negative  ion  densities  for  daytime 

and  nighttime.  The  positive  ion  densities  are  equal  to  the  sum  of 
the  electron  and  the  negative  ion  densities.  Curves  (a)  and  (b) 
are  the  calculated  altitude  profiles  for  electrons  for  day  and 
night,  respectively.  The  corresponding  day  and  night  profiles  for 
the  negative  ions  are  given  by  curves  (c)  and  (d).  Curves  (e)  and 
(f)  are  the  profiles  for  electron  and  negative  ions  calculate. 1  For 
nighttime  including  the  production  of  ions  from  energetic 
e lent  runs , 

Figure  8.  Calculated  altitude  profiles  of  negative  ions  for  day  (A)  and 

night  (B).  The  curves  are  O2  ( - ),  CO3  ( - ),  CO4  (***),  ND3 

(-•-*-),  NO2 —  (-0-0-),  0“  (-A-A-). 

Figure  9.  Calculated  altitude  profiles  of  02+  ("a”  for  day,  "b"  for  night) 
and  N0+  ("c"  for  day  and  "d"  for  night). 

Figure  10.  Calculated  day  (A)  and  night  (B)  altitude  profiles  of  the  dominant 

positive  ions  below  100  km.  The  curves  are  H  . (820)4  ( - ), 

N0+.(H20)  ( - ),  NO+.(H20)2  (•••),  N0+.(C02)  (-•  —  -),  H+(  H2(0  3 

(-0-0-). 

Figure  11.  Schematic  diagram  for  the  chemical  scheme  of  sodium  neutral  and 
ionic  species. 
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Figure  12.  Calculated  diurnal  variation  of  the  Na+  density  at  (a)  90  km,  (b) 
94  km,  and  (c)  100  km. 

Figure  13.  Calculated  altitude  profiles  for  Na  and  Na+  at  noon. 

Figure  14.  Calculated  altitude  profiles  of  atomic  sodium  using  different 

rates  for  reaction  (4-4).  The  curves  correspond  to  values  of  K* 

2  x  1 0  3 1  cm^  s  1  ( - ),  K*  =  7xl033  cm^  s  1  ( ),  K*  =*  4  x 

10  32  cm6  s  2  (-x-),  and  the  case  where  no  Na+  is  present  (-o-o-) 


Figure  15.  Calculated  profiles  of  the  Na  layer,  assuming  vertical  ion  drift 
velocities  of  0  ( - ),  -4  (-A-A-),  -10  (-o-o-),  and  -20  (-x-x-) 


Schematic  diagram  for  the  chemical  scheme  of  the  positive  ions  in 
the  D  and  lower  E  regions.  The  ions  are  grouped  into  two  families 
with  02+  and  N0+  being  the  primary  ions  in  each  group. 


69 


ALTITUDE  (KM)  ALTITUDE  (KM) 


105 

( e  )‘* 


DAY  SOURCES 

\  A 

c)  \  \(o) 

\t 

h 

Jr 


r — / 

•.  A 


NIGHT  SOURCES 


I  V 

(  c )!  \( b  ) 

I  \ 

I  \ 

j  j 

_ _ 

"  J 
/ 


°(f) 


DAY  SOURCES 


A'  If  - 


( a  )  A*' 


*  80 


..•*// 

//(c) 


LxJ 

§  75 

t- 


(  e )  /  // 


...•V/Md) 

•*  /  / 

/  / 

/  / 


10  10l  102  I03  10* 

PRODUCTION  (CM”3S~1 ) 


I0'2  10*1  10  101  10* 


PRODUCTION  (CM'3Sm) 


Figure  3.  The  production  rates  for 
positive  ions  for  daytime  and 
nighttime  conditions:  A)  Day 
(45°  Solar  Zenith  Angle  (SZA)), 
altitude  range  90-120  km);  B) 
Day  (45°  SZA),  altitude  range 
60-90  km;  c)  Night  (180°  SZA), 
altitude  range  60-120  km.  An 
altitude-independent  production 
rate  of  10  3  cm3  s_1  due  to 
gallactic  x-rays  is  also  in¬ 
cluded  in  the  model.  The 
curves  are:  a)  x-ray,  b) 
Lyman-g,  c)  Lyman-a ,  d)  ener¬ 
getic  electrons,  e)  ionization 
of  02(*A),  and  f)  cosmic  rays. 
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Figure  5.  Calculated  diurnal  variation  of  electron  density.  Profiles  are 
calculated  at  (a)  86  km,  (b)  90  km,  (c)  94  kra,  and  (d)  100  km. 
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Figure  7.  Calculated  total  electron  and  negative  ion  densities  for  daytime 

and  nighttime.  The  positive  ion  densities  are  equal  to  the  sum  of 
the  electron  and  the  negative  ion  densities.  Curves  (a)  and  (b) 
are  the  calculated  altitude  profiles  for  electrons  for  day  and 
night,  respectively.  The  corresponding  day  and  night  profiles  for 
the  negative  ions  are  given  by  curves  (c)  and  (d).  Curves  (e)  and 
(f)  are  the  profiles  for  electron  and  negative  ions  calculated  for 
nighttime  including  the  production  of  ions  from  energetic 
electrons . 
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Calculated  altitude  profiles  of  negative  ions  for  day  (A)  and 

night  (B).  The  curves  are  O2  ( - ),  CO3  ( - ),  COt+  (*’*), 

NO2  ~  (-0-0-),  0~  (-A-A-). 
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Figure  10-  Calculated  day  (A)  and  night  (B)  altitude  profiles  of  the  dominant 

positive  ions  below  100  km.  The  curves  are  H  .  (H20)i+  ( - ), 

N0+.(H20)  ( - ),  N0+. (H2 0)2  (***),  N0+.(C02)  ( - ),  H+(H20)3 


PART  IV 


Seasonal  Variation  of  the  Chemistry  and 
Transport  of  Sodium  Species  in  the  Mesosphere 


ABSTRACT 


Seasonal,  variations  have  been  observed  in  the  total  column  abundance  of 
Na  and  the  D-line  emission  rate.  Previous  studies  by  Garcia  and  Solomon 
(1985)  have  successfully  explained  seasonal  variations  of  0  and  O3  in  terms  of 
changes  in  the  vertical  diffusion  coefficient  K./z  induced  by  corresponding 
variations  in  the  breaking  of  gravity  waves  in  the  mesosphere.  We  have  incor¬ 
porated  transport  effects  into  our  model  and  investigated  the  sensitivity  of 
our  results  to  seasonal  changes  in  Kzz  and  background  atmospheric  density  and 
temperature.  Seasonal  changes  in  the  background  atmosphere  induce  a  winter 
maximum  in  the  Na  column  density  a  factor  of  1.4  higher  than  the  summer  mini¬ 
mum  in  qualitative  agreement  with  obser  vat  ions .  Seasonal  changes  in  Y^7?  in¬ 
duce  a  seasonal  variation  of  factors  of  3-4  in  the  Na  column  abundance  and  the 
D-line  emission,  peaking  at  the  equinoxes.  This  result  is  in  agreement  with 
low-latitude  D-line  observations,  but  in  disagreement  with  high-latitude  D- 
l ine  observations  at  all  latitudes.  We  cannot  thus  reconcile  the  Na  observa¬ 
tions  with  the  seasonal  changes  in  Kz?  within  the  context  of  the  present  pho¬ 
tochemical  model.  The  calculated  diurnal  varation  of  the  Na  column  density  is 
also  sensitive  to  seasonal  changes  in  background  trace  species,  particularly  0 
and  H2O. 
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INTRODUCTION 


The  calculations  presented  in  Part  I  concentrated  on  the  diurnal  varia¬ 
tion  of  background  and  sodium  species  in  the  mesosphere.  As  indicated  in  the 
discussion  of  the  available  sodium  data  (Part  I,  Section  2),  seasonal  patterns 
have  been  detected  in  both  the  Na  column  density  and  the  nighttime  sodium 
D-line  emission.  The  observed  Na  column  density  exhibits  a  maximum  in  winter 
at  all  latitudes  with  the  winter-summer  ratio  increasing  with  latitude 
(Simonich  et  al.,  1979;  Gibson  and  Sanford,  1971).  Spectral  analysis  of  the 
northern  hemisphere  data  on  D-line  emission  indicates  a  dominant  semiannual 
component  with  peaks  at  the  equinoxes  at  low  latitudes  while  the  high-latitude 
data  contains  a  stronger  annual  component  peaking  in  winter  (Fukuyama,  1977). 

Seasonal  variability  has  been  observed  or  deduced  for  other  trace  gases 
in  the  mesosphere.  Most  important,  the  O3  concentrations  deduced  from  the  IR 
observations  by  the  Solar  Mesosphere  Explorer  (SME)  exhibit  a  maximum  during 
early  fall  in  the  region  near  80  km  (Thomas  et  al.,  1984;  Garcia  and  Solomon, 
1986).  A  similar  maximum  is  observed  in  the  green  line  emission  of  atomic 
oxygen,  originating  near  100  km  (Garcia  and  Solomon,  1985). 

Assessment  of  the  mechanisms  proposed  to  explain  the  seasonal  variability 
of  0  and  O3  may  shed  some  light  on  the  processes  controlling  sodium  species. 

We  note  that  modeling  of  the  seasonal  behavior  of  atmospheric  species  may  re¬ 
quire  inclusion  of  meridional  transport  and  thus  demands  utilization  of  a  two- 
dimensional  model.  However,  recent  calculations  by  Garcia  and  Solomon  (1985) 
have  explained  the  observed  raidlatitude  seasonal  variations  in  0  and  O3  in 
terras  of  corresponding  changes  in  the  vertical  diffusion  coefficient  K zz . 

These  authors  have  introduced  a  parameterization  of  Kzz  based  on  the  breaking 
of  gravity  waves  in  the  mesosphere.  The  seasonal  variability  of  Kzz  is  in¬ 
duced  by  seasonal  changes  in  the  stratospheric  zonal  wind.  Although  their 
calculations  utilize  a  two-dimensional  model,  the  variability  in  0  and  O3  is 
controlled  primarily  by  vertical  diffusion.  These  results  suggest  an  initial 
analysis  of  the  seasonal  behavior  of  the  sodium  species  utilizing  a  1-D  model 
and  considering  seasonal  changes  in  the  eddy  coefficient,  background 
atmosphere,  and  rates  of  chemical  reactions. 

We  present  in  this  part  results  of  our  calculations  of  the  seasonal  vari¬ 
ation  of  background  and  sodium  trace  species  in  the  mesosphere  and  lower  ther¬ 
mosphere.  This  problem  requires  inclusion  of  a  self-consistent  calculation  of 
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"slow"  species  (i.e.,  species  affected  by  transport)  in  our  code.  The  calcu¬ 
lated  seasonal  behavior  of  the  mesosphere  has  also  implications  for  the  diur¬ 
nal  variation  of  the  Na  column  abundance  discussed  in  Part  I,  since  the  diur¬ 
nal  behavior  of  sodium  species  is  sensitive  not  only  to  the  adopted  chemical 
scheme  but  also  to  the  abundance  of  trace  gases  such  as  0  and  H2O. 

The  code  utilized  to  calculate  both  "fast"  and  "slow"  species  is  de¬ 
scribed  in  Section  2.  Section  3  presents  results  of  the  calculated  seasonal 
variation  for  background  trace  gases  (3(a))  and  sodium  species  (3(b)).  Sec¬ 
tion  4  addresses  again  the  diurnal  variation  of  the  sodium  column  density  ta¬ 
king  into  account  the  variability  of  the  background  atmosphere.  Conclusions 
are  presented  in  Section  5. 

2 .  TRANSPORT  CALCULATION 

Our  previous  calculations,  presented  in  Part  I,  have  fixed  the  concentra¬ 
tions  of  long-lived  species  to  the  values  of  Allen  et  al.  (1984).  The  effects 
of  eddy  diffusion  were  incorporated  in  our  diurnal  model  by  means  of  a  diffu¬ 
sion  term  consistent  with  the  assumed  profiles  for  slow  species.  Details  of 
this  procedure  were  given  in  Part  I  and  Rodriguez  et  al.  (1984).  This  approach 
is  satisfactory  for  studying  the  diurnal  behavior  of  short-lived  species  for 
fixed  conditions.  However,  it  precludes  the  ability  to  self-consistently  cal¬ 
culate  concentrations  for  arbitrary  vertical  eddy  diffusion  coefficient  (Kzz) 
and  solar  zenith  angle.  The  need  for  such  calculations  is  particularly  im¬ 
portant  in  studying  the  effect  of  changes  in  Kzz  on  the  column  abundance  of 
sodium. 

We  extended  our  one-dimensional  mesospheric  model  to  include  the  effects 
of  eddy  and  molecular  diffusion.  The  species  calculated  by  the  model  are: 
odd-oxygen  species,  0X  —  0,  0(1D),  O3;  odd-hydrogen  species,  H0X  —  H,  OH, 

HO2 ,  H2O2;  sodium  species,  NaX  —  Na,  NaO,  Na02,  NaOH,  and  H2O,  H2,  and  0  2* 

The  nitrogen  species  are  not  included  in  this  calculation,  since  their  effect 
on  0  ,  H0k,  and  NaX  is  minor.  The  photochemical  lifetimes  of  the  odd-oxygen 
and  hydrogen  species  are  less  than  one  cay  at  altitudes  below  85  km.  However, 
these  lifetimes  increase  rapidly  above  85  km  and  the  concentrations  of  0X  and 
H0X  are  then  primarily  determined  by  transport  processes.  The  principal 
sources  of  H0X  and  0X  are  photo'iysis  of  H2O  and  02>  respectively.  A  large 
fraction  of  the  H0X  produced  by  H2O  photolysis  above  75  km  recombines  to  form 
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H2  (Liu  and  Donahue,  1974).  We  thus  have  to  calculate  all  the  above  species 
in  a  formulation  including  eddy  and  molecular  transport. 

The  model  solves  the  time-averaged  continuity  equation  for  the  flux  4^  of 
the  ic^  species 


d<4>i> 


=  <?L>  ~  <Lt>, 


where  and  L^  represent  the  production  and  loss  processes  for  the  i11  spe¬ 
cies.  The  brackets  <  >  denote  an  average  over  a  full  diurnal  cycle.  Equation 
(1)  translates  into  an  equation  for  the  densities  n^  by  using  the  expression 
for  the  vertical  flux  (Hunten,  1975) 


<<ti>  =  -Kn 


d<fi> 


-  Di  [ 


l  d<nj>  mjg  (1  +  aj)  dx1 
n  dz  kT  T  dz  - 


In  the  above,  K  and  D^  denote  the  eddy  and  molecular  diffusion  coefficients, 
respectively;  n,  the  background  atmospheric  density;  ,  the  molecular  mass  of 
species  i;  and  a ^ ,  the  thermal  diffusion  coefficient.  Equations  (1)  and  (2) 
can  be  combined  into  a  single  second-order  differential  equation  for  the  spe¬ 
cies  concentration  n^  with  unique  solutions  for  given  boundary  conditions  at 
the  bottom  and  top  altitudes  (60  and  120  km). 

The  solution  of  equations  (1)  and  (2)  requires  specifying  the  production 
and  loss  mechanisms,  the  boundary  conditions,  and  the  eddy  and  molecular  dif¬ 
fusion  coefficients.  Each  of  these  will  be  discussed  in  Sections  2(a)  through 
2(c).  The  actual  computational  procedure  for  the  solution  is  presented  in 
Section  2(d). 

( a )  Production  and  Loss  Mechanisms 

The  production  and  loss  mechanisms  for  the  long-lived  species  are 
obtained  from  the  chemical  scheme  presented  in  our  diurnal  model  (Part  I)  and 
summarized  in  Table  1.  In  the  lower  mesosphere,  H0X  is  produced  mainly  by 

H2 0  +  0(!d)  -*■  2 OH  (3) 

and  removed  by 

OH  +  HO2  +  H2O  +  O2 


OH  +  OH  +  H2O  +  0. 


0KW5 


The  production  of  H0X  above  70  km  occurs  primarily  through  photolysis  of 


water 

H20  +  hv  -*•  OH  +  H.  (5) 

An  additional  sink  of  H0X  in  this  region  is  the  reaction 

H  +  H02  H2  +  O2  (6a) 

+  H2  0  +  0 .  (60) 

Reactions  (5)  and  (6a)  eventually  convert  all  water  into  H0x  (mostly  H)  and 
H2 .  Molecular  hydrogen  is  the  dominant  hydrogen-bearing  species  below  120  km. 

Odd  oxygen  is  produced  by  photolysis  of  02: 

02  +  hv  20;  (7) 

the  main  sinks  for  0x  are  the  reactions 

0  +  OH  *  02  +  H 
H  +  O3  +  OH  +  O2 
0  +  H02  +  OH  +  02  . 

Above  90  km,  three-body  recombination  of  0  becomes  important: 

O+O+M+O2+M.  (11) 

Reactions  (7)  -  (11)  also  constitute  the  main  sinks  and  sources  of  02* 

We  assume  that  the  sodium  species  NaX  are  produced  by  meteor  ablation 
between  80  and  100  km  (Gadsden,  lq68;  Hunten  et  al.,  1^80).  These  species 
diffuse  down  to  the  troposphere  where  they  are  removed  by  washout  with  a  time 
constant  of  10  days  (Liu  and  Reid,  1979).  We  have  adopted  the  production 
profile  of  Hunten  et  al.  (1980),  shown  in  Figure  1. 

(b)  Boundary  Conditions 

Equations  (1)  and  (2)  are  solved  between  60  and  120  km.  We  adopt 
a  fixed  mixing  ratio  boundary  condition  for  ali  jpecies  at  60  km.  The  mixing 
ratios  of  H20  (7  ppmv),  H2  (0.5  ppmv),  and  02  (21%)  at  60  km  are  <.dopted  from 
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previous  stratospheric  calculations.  At  this  altitude,  both  0x  and  H0X  have 
photochemical  lifetimes  of  less  than  one  day.  The  diurnal  average  of  their 
mixing  ratio  is  then  independent  of  transport  and  can  be  obtained  from  our 
diurnal  calculation. 

We  adopt  flux  boundary  conditions  at  the  upper  boundary.  Fluxes  are 
assumed  to  be  zero  for  all  slow  species  except  as  indicated  below.  For  the 
oxygen  species, 

2<f>  (120  km)  =  (120  km)  =  /“  2J7  [02J  dz,  (12) 

2  x  120  km 

where  J7  is  the  photolysis  rate  in  reaction  (7).  Equation  (2)  neglects 
recombination  of  0  above  120  km.  The  integral  in  equation  (13)  can  be 
approximated  by 


/"  J7  [0,  1  dz  «*  J7  ( 120  km)  n  (120  km)  H  (120  km),  (13) 

120  km  °2  °2 

where  Hq2  is  the  scale  height  of  O2  (approximately  6.5  km). 

Atomic  hydrogen  escapes  from  Earth  at  a  rate  of  2-3  x  108  cm'2  s_1  (Hun- 
ten  and  Donahue,  1976).  We  therefore  have,  at  120  km, 

+  2(f)  +  2<t>  =  2-3  x  108  cm-2  s_1.  (14a) 

nu^  ri2  U  1^2 

Previous  models  by  Liu  and  Donahue  (1974)  show  that  most  of  the  elemental  hy¬ 
drogen  flux  at  120  km  is  in  the  form  of  H2.  We  thus  adopt  the  boundary  condi- 
t  ion 


4>H0  (120  km)  -  $h?o(  120  km)  =0 

x  z  ( 14b) 

(120  km)  =  1-1.5  x  108  cm-2  s— 1 . 

(c)  Molecular  and  Eddy  Diffusion  Coefficients 

Turbulent  mixing  in  the  vertical  direction  is  parameterized  by  an 
eddy  diffusion  coefficient  Kzz  (Hunten,  1975).  Several  authors  have  proposed 
different  values  of  Kzz  to  explain  the  observed  altitude  distribution  of  meso¬ 
spheric  species  (Hunten,  1975;  Allen  et  al.,  1981;  Garcia  and  Solomon,  1985). 
While  values  such  as  those  of  Allen  et  al.  (1981)  are  adjusted  to  fit  specific 
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observations,  the  parameterization  of  Garcia  and  Solomon  (1986)  relates 
turbulent  mixing  to  physical  processes  in  the  mesosphere  (i.e.,  gravity-wave 
breaking).  We  are  thus  particularly  interested  in  the  variations  induced  by 
the  seasonally-dependent  Kzz  proposed  by  Garcia  and  Solomon  (1985).  Profiles 
for  these  eddy  diffusion  coefficients  are  shown  in  Figure  2. 

The  effects  of  molecular  diffusion  become  important  above  100  km.  Molec¬ 
ular  diffusion  coefficients  are  given  by  (Banks  and  Kockarts,  1973) 

b. 

Di  =  —  (15a) 

n 

where 

=  Al  TSi,  (15b) 

and  n  and  T  are  the  background  atmospheric  density  and  temperature.  Table  2 
presents  our  adopted  values  for  the  parameters  A^  and  s^  for  each  species. 
Adopted  values  for  the  thermal  diffusion  coefficient  are  -0.27,  -0.39,  and 
-0.31  for  0X,  H0X,  and  K2 ,  and  zero  otherwise  (Allen  et  al.,  1981). 

Total  densities  and  temperature  of  the  background  atmosphere  exhibit 
considerable  seasonal  variability  between  80  and  100  km  (see,  i.e.,  Cole  and 
Kantor,  1978).  We  have  adopted  the  seasonally-dependent  model  atmosphere  of 
Cole  and  Kantor  (1978)  (z  <  90  km)  and  CIRA,  1972)  (z  >  90  km)  in  the 
calculations  described  below,  unless  otherwise  specified.  The  seasonally- 
dependent  background  temperature  and  density  adopted  in  our  calculations  are 
shown  in  Figures  3  and  4. 

(d)  Computational  Procedures 

The  basic  structure  of  the  coupled  diurnal-diffusion  code  is 
illustrated  by  the  flow  chart  in  Figure  5.  An  initial  guess  for  the  altitude 
profile  of  slow  species  is  given  as  input.  Our  previous  diurnal  code  is  then 
used  to  solve  for  the  diurnal  variation  of  fast  species  consistent  with  the 
fixed  slow  species  density  at  each  altitude.  The  calculated  diurnal  species 
are  then  used  to  calculate  the  time-averaged  <P^>  and  <L^>  for  equations  (l) 
and  (2).  Solutions  of  these  equations  give  us  a  new  estimate  for  the  slow 
species  which  is  used  again  as  Lnput  to  the  diurnal  code.  The  process  is 
iterated  until  the  difference  in  the  densities  calculated  in  two  successive 
iterations  is  less  than  some  specified  value. 
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With  the  assumed  values  for  Kzz  and  1)^,  resident  times  due  to  eddy  and 
molecuLar  diffusion  are  of  the  order  of  days,  except  at  night  between  80  and 
85  km  (Allen  et  al.,  1984).  In  this  region,  the  effect  of  transport  is 
introduced  by  means  of  a  constant  term  P*  in  the  diurnal  continuity  equation 
(Part  I  and  Rodriguez  et  al.,  1984). 


3.  SEASONAL  VARIATIONS 

The  coupled  diurnal  transport  model  described  in  the  previous  section  was 

used  to  calculate  seasonally-dependent  densities  of  the  background  trice  gases 

and  sodium  species.  Calculations  were  performed  for  spring,  summer,  fall,  and 

winter  conditions  at  45°  latitude  using  the  respective  Kzz  profiles  and 

background  atmospheres.  Results  are  presented  for  fall  and  winter  to 

illustrate  the  effect  of  maximum  and  minimum  values  of  K„  . 

zz 

Our  results  for  background  trace  species  are  presented  in  Section  3(a). 
Section  3(b)  discusses  our  calculations  for  the  seasonal  variation  of  the 
sodium  column  abundance  and  n-line  emission. 

(a)  Background  Trace  Species 

Results  of  our  calculations  for  H2O,  0  (noon),  and  H  (noon)  are 
presented  in  Figures  4  through  6  for  fall  and  winter  conditions.  These  fall 
and  winter  trace  species  were  calculated  by  allowing  for  changes  in  the  photo¬ 
lysis  rates,  background  atmosphere  (n,T),  and  eddy  diffusion  coefficients.  We 
also  include  for  comparison  the  values  used  in  the  calculations  described  in 
Parts  1  and  ill,  corresponding  to  the  profiles  calcualted  by  of  Allen  et  al. 
(1984)  which  used  an  adjusted  Kz?  to  simulate  the  observed  concentrations. 

Water  is  photodissociated  primarily  by  Lyman-alpha  radiation  above 
75  km.  The  abundance  of  water  in  the  upper  mesosphere  is  then  determined  by 
the  efficiency  of  the  vertical  transport  from  the  lower  atmosphere.  The  eddy 
diffusion  coefficients  derived  by  Garcia  and  Solomon  (1985)  are  considerably 
lower  in  fall  than  in  winter.  As  a  result,  water  is  transported  less 
efficiently  to  the  upper  mesosphere  in  the  fall,  leading  to  considerably  lower 
H2O  abundances.  This  seasonal  betiavior  is  illustrated  in  Figure  6. 

The  H  densities  shown  in  Figure  7  also  exhibit  lower  values  in  the  tall, 
particularly  above  85  km.  In  this  region,  most  of  the  MO  species  are  in  the 
form  of  atomic  hydrogen.  Since  the  source  of  H0x  is  photolysis  of  H ,  the 


the  1<  r  water  concentration  in  the  fall  leads  to  lower  H0x  densities. 

Higher  photolysis  rates  for  H2O  in  fall  tend  to  partially  compensate  for  the 
lower  H2O  concentration  resulting  In  '.iOx  density  differences  of  less  than  50 
percent.  The  variation  below  85  vra  reflects  the  effect  of  changes  in  solar 
zenith  angle  resulting  in  larger  preduction  rate  and  higher  concentration  of 
hydrogen  in  the  fall. 

Figure  8  shows  calculated  0  densities  which  exhibit  larger  values  in  fall 
than  in  winter  above  85  km.  This  behavior  is  a  consequence  of:  1)  the  lower 
concentrations  of  H0X  species  which  constitute  a  sink  of  0X  through  the 
reactions  of  OH  and  H  with  03,  and  2)  the  slower  transport  of  atomic  oxygen, 
produced  by  photolysis  of  O2  above  90  km,  out  of  this  region.  Calculated 
densities  of  O3  at  midnight  are  shown  in  Figure  9  for  fall  and  winter.  The  O3 
profile  exhibits  the  minimum  near  80  km  typical  of  both  observations  and  cal¬ 
culations  (see  Allen  et  al.,  1984).  The  seasonal  changes  in  0,  H,  and  O3 
illustrated  in  Figures  5  through  9  are  consistent  with  the  midlatitude  2-D 
results  of  Garcia  and  Solomon  (1985). 

( b )  Seasonal  Variation  of  Sodium  Abundance  and  D-llne  Emission 

Observations  of  the  mesospheric  sodium  layer  have  established  the 
existence  of  a  maximum  in  the  Na  column  density  during  winter.  This  maximum 
is  a  factor  of  two  larger  than  the  summer  minimum  at  low  latitudes  (Simonich 
et  al.,  1979)  and  up  to  three  to  four  times  larger  at  high  latitudes  (Gibson 
and  Sandford,  1971).  Three  different  processes  could  contribute  to  the  ob¬ 
served  seasonal  variation  of  the  sodium  column  density:  changes  in  the  ver¬ 
tical  transport,  changes  in  photochemical  control  (Swider,  1985),  and  changes 
in  downward  fluxes  of  Na+  from  above  95  Km.  The  relative  contributions  of 
each  of  the  above  to  the  seasonal  variation  have  been  assessed  with  our  one- 
dimensional  model. 

We  can  best  understand  the  direct  effect  of  vertical  transport  by 
considering  the  behavior  of  a  species  with  a  source  restricted  to  the  upper 
atmosphere  (z  >  z\  )  and  a  sink  restricted  to  the  lower  atmosphere  (z  <  Z2) 
with  z\  >  Z2 .  The  mixing  ratio  f  at  an  altitude  z  >  Z2  is  then  given  by 


f  (z)  =  fr 


/ 


z _ 4>(z’  ) 

z>  K77 (z ' )n( z  '  ) 


dz  ' 


(16) 


where  fo  is  the  mixing  ratio  at  Z2 ,  Kzz(z)  is  the  vertical  eddy  diffusion 


coefficient,  n  is  the  background  atmospheric  density,  and  <f>(z)  is  the  downward 
flux  (cm-2  s'1)  equal  to  the  integrated  production  J*  Pdz  in  the  upper  atmo¬ 
sphere.  Equation  (16)  can  be  applied  to  the  family  of  sodium  species  NaX  = 

Na  +  NaO  +  NaC>2  +  NaOH,  which  is  produced  by  meteor  ablation  between  80  and 
100  km  (Gadsden,  1968;  Hunten,  1981)  and  lost  in  the  stratosphere  or  tropo¬ 
sphere  by  attachment  to  aerosol  particles,  polymer  formation,  or  washout 
(Hunten,  1981;  Lamb  and  Benson,  1986;  Liu  and  Reid,  1979).  We  can  also  dis¬ 
regard  fo  in  (16)  for  Z2  <  30  km  when  compared  to  the  integral  term  for 
altitudes  below  the  Na  peak  at  90  km.  We  can  then  see  that  the  abundance  of 
NaX  (and  thus  of  Na)  in  the  sodium  layer  (80  to  100  km)  will  be  inversely 
proportional  to  the  adopted  eddy  diffusion  coefficient.  The  largest  contri¬ 
butions  to  the  integral  in  (16)  come  from  altitudes  near  the  region  of  meteor 
ablation,  thus  it  is  most  sensitive  to  Kzz  between  80  and  100  km.  Considera¬ 
tion  of  the  eddy  diffusion  coefficients  in  Figure  2  in  conjunction  with  (16) 
suggests  that  the  mixing  ratio  of  NaX  should  have  maxima  at  equinoxes  with 
minima  at  the  solstices. 

Another  important  seasonal  effect  has  been  pointed  out  by  Swider  (1985). 
The  density  of  Na  in  the  bottomside  layer  is  controlled  by  the  rapid  formation 
of  Na02  through 

Na  +  02  +  M  Na02  +  M.  (17) 

The  seasonal  dependence  of  the  temperature  and  background  density  below  90  km 
combine  to  minimize  the  loss  of  Na  through  (17)  in  the  winter,  thus  maximizing 
the  sodium  density  in  this  season.  We  have  investigated  the  above  effect  by 
incorporating  the  seasonally-dependent  background  atmosphere  of  Cole  and 
Kantor  (1978)  into  our  one-dimensional  model. 

The  column  density  of  Na  could  also  be  sensitive  to  changes  in  the 
topside  scale  height  due  to  seasonal  variations  of  Na+.  Hanson  and  Donaldson 
(1967)  suggested  a  mechanism  to  explain  the  small  topside  scale  height,  by 
which  upward  diffusion  of  neutral  Na  should  be  roughly  balanced  by  the 
downward  ion  drift  flux  W[Na+],  where  W  is  the  ion  drift  velocity.  We  can 
then  expect  that  changes  in  either  W  or  [Na+]  will  change  the  upward  flux  of 
Na,  and  thus  the  topside  scale  height  and  total  column  density.  Details  of 
this  effect  were  discussed  in  Part  III.  We  have  assumed  a  constant  vertical 


drift  velocity  of  -20  cm/s  as  in  Part  III  and  investigated  the  sensitivity  of 
our  results  to  different  diurnally-averaged  Na+  resulting  from  changing  length 
of  day. 

The  sensitivity  of  the  calculated  Na  column  abundance  to  each  of  the 
above  parameters  is  illustrated  by  the  results  in  Table  3.  The  one-dimension¬ 
al  model  was  run  for  conditions  typical  of  45°N  spring,  summer,  fall,  and  win¬ 
ter.  In  the  first  set  of  runs  (Case  1),  the  eddy  diffusion  coefficient  was 
held  constant  at  the  summer  values  in  Figure  2,  while  the  temperature  and 
background  density  were  varied  according  to  Cole  and  Kantor  (1978)  and  CIRA 
(1972).  The  second  set  of  runs  (Case  2)  held  the  background  atmosphere 
constant  while  varying  Kzz  according  to  Figure  2.  The  third  set  (Case  3) 
changed  both  the  eddy  diffusion  coefficient  and  the  background  atmosphere. 

For  each  set,  we  also  investigated  the  effect  of  excluding  and  including  the 
vertical  ion  drift. 

The  results  in  Table  3  indicate  that  a  fall-winter  maximum  is  obtained  if 
we  vary  the  background  atmosphere  (Case  1).  The  winter-summer  contrast  (~ 

1.4)  is,  however,  closer  to  the  low-latitude  value  of  ~  2.0  of  Simonich  et 
al.,  1979  (23°  S)  than  to  the  high  latitude  value  of  3  to  4  observed  by  Gibson 
and  Sandford  (1971)  at  51°  N.  We  can  expect  the  seasonal  effect  in  Case  1  to 
increase  with  latitude,  in  response  to  a  larger  seasonal  variation  of  the 
background  atmosphere,  and  larger  contrast  in  the  length  of  day  between  summer 
and  winter. 


The  results  in  Cases  2  and  3  indicate  that  the  response  to  seasonal 
variations  of  Kzz  is  much  larger  than  the  corresponding  response  to  variations 
in  the  background  atmosphere.  The  column  density  of  Na  changes  by  about  a 
factor  of  3  between  minimum  and  maximum.  We  note,  however,  that  the  results 
exhibit  two  maxima  at  the  equinoxes,  whereas  observations  indicate  a  single 
maximum  during  winter. 


Our  results  indicate  that  the  observed  winter  maximum  in  Na  column 
density  can  be  obtained  from  changes  in  tiie  background  temperature  and  density 
(Swider,  1983).  Although  inclusion  of  ions  changes  the  absolute  magnitude  of 
the  column  density,  it  does  not  seem  to  affect  the  seasonal  contrast.  The 
magnitude  of  the  variation  is,  however,  smaller  than  observed.  The  proposed 
behavior  of  Kzz  by  Garcia  and  Solomon  (1985)  thus  seems  inconsistent  with  the 
observed  seasonal  variation  of  Na ,  at  least  within  the  context  of  a  1-D 
model.  Further  work  is  needed  to  reconcile  the  O3,  0,  and  Na  behavior. 


In  contrast  to  the  lidar  measurements,  the  Na  D-line  emission  observa¬ 
tions  do  not  give  direct  measurement  of  the  column  abundance.  The  D-line 
emission  is  given  by: 

Iq(R)  =  /  k2[Na0][0]  -  /  kj  lNaJ[03],  (18) 

where  k^  and  k2  are  the  reaction  rate  constants  for  the  reactions  Na  +  0-j  + 

Na  0  +  0  and  Na  0+0  -*■  Na  +  02,  respectively  (Swider,  1986).  It  should  be 
noted  that  the  magnitude  of  the  intensity  of  the  sodium  D  line  depends  on  the 
convolution  of  the  sodium  and  ozone  altitude  profiles.  Since  O3  decreases 
rapidly  with  height  above  80-85  km,  the  D-line  emission  is  most  sensitive  to 
the  Na  concentration  in  the  bottomside  layer. 

The  seasonal  behavior  of  the  observed  intensity  of  the  sodium  night-glow 
is  shown  in  Fig.  10.  Emission  rates  in  this  figure  are  normalized  to  50  R. 
This  emission  exhibits  maxima  at  the  equinoxes  and  minima  at  the  solstices 
(Kirchoff  et  aL.,  1981)  at  a  latitude  of  23°S.  For  high-latitudes,  however, 
the  D-line  nightglow  exhibits  a  maximum  in  winter,  with  lesser  peaks  at  the 
equinoxes  (Fukuyama,  1977).  This  behavior  of  the  sodium  D  line  cannot  be 
explained  solely  by  seasonal  changes  in  the  J  rates  as  these  changes  would 
tend  to  produce  a  maximum  in  summer.  However,  the  low-latitude  behavior  is 
consistent  with  the  expected  response  of  Na  to  the  effect  of  eddy  transport 
where  the  larger  values  of  K  during  solstices  provides  more  ef f icient 
removal  of  Na  by  transport  into  the  stratosphere. 

Our  calculated  emission  rates  for  each  season  are  also  shown  in  Figure 
10.  The  absolute  magnitude  of  the  emission  rate  is  proportional  to  the  total 
NaX  production  and  the  branching  for  Na(zF)  production  by  (18);  since  these 
quantities  are  uncertain,  we  concentrate  on  the  seasonal  changes  rather  than 
the  absolute  magnitude  (histogram).  The  seasonal  behavior  shown  in  Fig.  10  is 
In  good  agreement  with  observation  at  23°S.  The  calculated  absolute  magnitude 
of  the  seasonal  change  is  larger  than  observed.  Our  calculations  do  not 
reproduce,  however,  the  winter  peak  observed  at  high  latitudes  (Fukuyama, 
1977). 

We  note  that  our  calculated  equinox  maximum  in  the  D-line  emission  is  a 
result  of  maxima  in  both  NaX  and  O3  .  A  smaller  equinox  maximum  would  still  be 
obtained  if  we  only  included  the  O3  density  behavior,  as  pointed  out  by  Swider 


(1986).  This  can  explain  the  D-line  observations  at  low  latitudes.  On  the 
other  hand,  the  winter  maximum  at  high  latitudes  is  probably  related  directly 
to  the  maximum  in  the  Na  abundance. 

4 .  CONCLUSIONS 

The  results  presented  above  have  utilized  a  1-D  model  and  relied  on  a 
particular  parameterization  of  Kzz  to  obtain  most  seasonal  effects.  Two- 
dimensional  effects,  however,  could  be  important.  The  values  of  Kzz  deduced 
by  Garcia  and  Solomon  (1985)  assume  a  particular  spectrum  for  breaking  gravity 
waves,  for  which  there  is  little  observational  data  at  present.  The  above 
uncertainties  in  our  calculation  imply  that  we  sould  concentrate  on 
qualititat ive  rather  than  quantitative  features  of  the  seasonal  variation  of 
different  species  in  the  evaluation  of  the  model  results. 

Our  model  calculates  seasonal  profiles  for  0x>  H0X,  and  H2O,  which  are 

consistent  with  the  results  of  Garcia  and  Solomon  (1985)  and  with  available 

observational  data.  The  seasonal  parameterization  of  K  accounts  for  the 

z  z 

first-order  features  of  the  seasonal  behavior  of  0X,  H0X,  and  H2O  in  the 
mesosphere . 

The  Kzz  parameterization  is  less  successful  in  describing  the  observed 
seasonal  behavior  of  the  Na  column  abundance.  Although  a  large  seasonal 
variation  is  obtained,  the  calculated  maxima  occur  at  the  equinoxes,  while 
observations  indicate  a  peak  in  winter.  The  seasonal  behavior  of  Na  is  better 
explained  by  the  variations  in  background  temperature  and  density  (Swider, 
1985)  if  we  disregard  changes  in  Kzz«  The  inconsistency  between  the  model 
required  to  explain  the  seasonal  behavior  of  0x  and  that  for  Na  is  not 
resolved  at  present.  Further  work  on  this  issue  would  provide  important 
insights  on  mesospheric  processes. 

Similar  conclusions  hold  for  the  results  on  the  Na  D-line  emission.  The 
agreement  between  our  calculations  and  the  low-latitude  observations  of 
Kirchoff  et  al.  (1981)  may  be  misleading.  The  same  qualitative  agreement 
could  be  obtained  by  involving  seasonal  changes  in  03.  Swider  (1986) 
indicated  that  most  of  the  D-line  emission  is  emitted  within  2  to  3  km  below 
the  Na  peak  at  90  km,  where  the  seasonal  variation  of  Na  is  expected  to  be 
small,  particularly  at  low  latitudes.  The  behavior  of  the  D-line  emission 
thus  decouples  from  that  of  the  Na  abundance.  At  the  same  time,  the  high- 
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latitude  seasonal  variability  exhibits  a  winter  maximum  (Fukuyama,  1977). 

This  may  be  due  to  a  stronger  seasonal  variation  of  Na  near  the  peak  at  these 
latitudes,  in  which  case  the  behavior  of  the  D-line  emission  and  the  Na 
abundance  are  coupled. 

An  extensive  data  base  exists  for  Na  density  and  emission  rates.  Since 
these  observations  are  ground-based,  they  can  be  carried  out  over  extensive 
periods  of  time  at  a  wide  range  of  locations..  The  sodium  observations  could 
thus  provide  an  important  complement  to  the  species  most  often  utilized  at 
present  in  validating  mesospheric  models  (i.e.,  O3,  0,  OH).  The  question 
raised  in  this  section  then  points  toward  important  areas  of  future  research. 

The  sensitivity  of  the  seasonal  variation  of  Na  to  changes  in  Kzz  is 
larger  than  that  of  O3  and  0.  It  is  thus  puzzling  that  the  present  transport 
parameterization  does  not  yield  the  observed  sodium  seasonal  behavior.  A  more 
realistic  study  requires  utilization  of  a  2-D  model  for  the  mesosphere. 
Although  beyond  the  scope  of  the  present  research,  such  study  can  provide 
additional  light  on  mesospheric  transport  mechanisms  of  importance  not  only 
for  sodium  species  but  also  for  H2O,  0X,  and  H0X.  The  parameterization  of  Kz? 
should  be  examined  more  closely  to  determine  possible  modifications  which 
could  explain  all  the  available  observations.  Finally,  chemical  mechanisms 
for  Na  species  not  included  in  this  study  should  also  be  considered  (see,  for 
example,  Hunten,  1981). 

Expansion  of  reanalysis  of  the  Na  data  base  could  also  shed  light  on  the 
issues  raised  above.  The  chemical  mechanism  proposed  by  Swider  (1985)  implies 
little  seasonal  variation  of  Na  near  the  peak,  whereas  changes  in  Kzz  would 
induce  changes  in  Na  at  all  altitudes.  Further  examination  of  data  for 
possible  detection  of  seasonal  changes  in  Na  peak  altitude  and  magnitude  is 
necessary.  Data  on  the  diurnal  variation  of  the  Na  column  abundance  should 
also  be  analyzed  for  possible  seasonal  trends.  The  extent  of  the  correlation 
between  Na  D-line  emission  and  Na  column  density  at  different  altitudes  should 
also  be  determined. 
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Table  1 

Production  and  Loss  Mechanisas  for  Diffusive  Species 


Species  Production 


H20  +  (X1!))  >  20H 

H2O  +  hv  +  H  +  OH 

OH  +  H02  +  H2O  +  02 

011  +  OH  +  H20  +  0 

H  +  H02 

H  +  H02 

+  H2  +  02 

-»■  H20  +  0 

02  +  hv  +  0  +  0 

H  +  O3  -*■  OH  +  O2 

0  +  H02  >  OH  +  02 

O  +  O  +  M+O2  +  M 

0  +  OH  ->■ 

O2  +  H 

OH  +  H02  +  H2  0  +  02 

OH  +  OH  ->•  H2O  +  0 

H20  +  OOD)  >  20H 

H20  +  hv  ■+■  H  +  OH 

H  +  H02  ■+  H2  +  02 

0  +  OH  -»  02  +  H 

H  +  03  +  OH  +  02 

0  +  ho2  +  oh  +  o2 

0  +  0  +  M  +  02  +  M 

O2  +  hv 

-*■  20 

Meteor  ablation 

Washout 

in  troposphere 

Table  2 

Molecular  Diffusion  Coefficients 


D  = 

b  s 
n 

ATS  2  _x 

=  -  craz  s  1 

n 

A 

s 

Reference 

H2 

2.6 

X 

1017 

0.750 

Mason  and  Marrero  (1970) 

H0x  (H) 

00 

• 

X 

1017 

0.7 

Banks  and  Kockarts  (1973) 

H2O 

1.37 

X 

1016 

1.07 

Mason  and  Marrero  (1970) 

O2 

6.84 

X 

1016 

0.77 

Estimated3 

\  (0) 

9.68 

X 

1016 

0.77 

Mason  and  Marrero  (1970) 

aEstimated  from  the  0  diffusion  coefficient  assuming  an  inverse  dependence  on 
the  square  root  of  the  molecular  mass  (Mason  and  Marrero,  1970). 


Table  1 

Calculated  Coluan  Densities  (ca-2)  of  Na  at  45°  N 


No  Na  + 


With  Na 
lW  =  -20  era/ s) 


Case 

1:  Constant 

Kzz* 

seasonal 

T. 

n 

Spring 

3.7 

X 

109 

2.8 

X 

10  9 

Summer 

3.1 

X 

109 

2.3 

X 

109 

Fall 

4.3 

X 

109 

3.3 

X 

109 

Winter 

4.0 

X 

109 

3.2 

X 

109 

Case 

2:  Seasonal 

Kzz. 

constant 

T, 

n 

Spring 

10.7 

X 

109 

8.1 

X 

109 

Summer 

3.9 

X 

109 

3.0 

X 

109 

Fall 

9.8 

X 

109 

7.4 

X 

109 

Winter 

2.7 

X 

109 

2.1 

X 

109 

Case 

3:  Seasonal 

Kzz  > 

T,  n 

Spring 

10.3 

X 

109 

7.6 

X 

109 

Summer 

3.1 

X 

109 

2.3 

X 

109 

Fall 

11.8 

X 

109 

8.9 

X 

109 

Winter 

3.3 

X 

109 

2.8 

X 

109 
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Captions 

1.  Adopted  production  of  Na  atoms  by  meteor  ablation.  The  shape  of 
the  profile  is  taken  from  Hunten  et  al.  (1980)  and  scaled  to  yield 
an  NaX  density  of  3  x  103  cm-3  at  90  km  for  the  summer  K  in 
Figure  2. 

2.  Vertical  eddy  diffusion  coefficient  K  adopted  in  our  calculations 

for  spring  ( - ),  summer  ( - ),  fall  (-•-),  and  winter  (•••) 

conditions.  These  profiles  are  taken  from  Garcia  and  Solomon's 
(1985)  parameterization  of  mixing  induced  by  gravity-wave  breaking 
at  45°  latitude.  The  profile  deduced  by  Allen  et  al.  (1981)  is 
also  shown  for  comparison. 

3.  Background  neutral  densities  adopted  in  our  calculations  for  spring 

( - ),  summer  ( - ),  fall  (***),  and  winter  (-•-)  conditions  below 

90  km  from  Cole  and  Kantor  (1978). 

4.  Same  as  Figure  3  but  for  the  background  temperature. 

5.  Flow  chart  for  coupled  diurnal-diffusion  program.  The  box  marked 
with  an  asterisk  (*)  corresponds  to  our  previous  diurnal  code.  The 
term  P£  represents  the  effects  of  transport  in  the  diurnal 
calculation,  as  explained  in  Rodriguez  et  al.  (1984). 

6.  Mixing  ratio  of  H2O  for  different  seasonal  conditions  used  in 
calculations.  The  seasonally-dependent  eddy  diffusion  coefficients 
of  Garcia  and  Solomon  (1985)  have  been  used  to  calculate  the  H2O 

mixing  ratio  for  fall  ( - )  and  winter  ( - )  conditions.  The 

profiLe  calculated  by  Allen  et  al.  (1984)  (-0-),  used  in  our  Fart  I 
calculations,  is  also  shown  for  comparison. 

7.  Same  as  Figure  4  but  for  the  density  of  H  at  noon. 

8.  Same  as  Figure  4  but  for  the  density  of  0  at  noon. 

9.  Density  of  O3  at  midnight,  calculated  For  fall  ( - )  and  winter 

( - )  conditions  using  the  K  parameter  t /.al  ton  >f  Garcia  and 

Solomon  (1985). 

10.  Calculated  seasonal  variation  of  the  sodium  D  line  nightglow  at 
45°N  (histogram)  compared  to  the  observed  seasonal  variation  at 

23°S  ( - )  (Kirchoff  et  al.,  1981)  and  4l°N  ( - )  (Fukuyama, 

1977).  Calculated  and  observed  emissions  have  been  normalized  to 

50R.  The  northern  hemisphere  data  of  Fukuyama  (1977)  ( - )  has 

been  shifted  by  six  months  for  comparison  with  the  southern 
seasons . 
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Figure  1.  Adopted  production  of  Na  atoms  by  meteor  ablation.  The  shape  of 

the  profile  is  taken  from  Hunten  et  al.  (1980)  and  scaled  to  yield 
an  NaX  density  of  3  *  10^  cm  3  at  90  km  for  th  esummer  K  in 
Figure  2.  zz 


Figure  2. 
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Vertical  eddy  diffusion  coefficient  K  adopted  in  our 

calculations  for  spring  ( - ),  summer  ( - ),  fall  (•••),  and 

winter  (-•-)  conditions.  These  profiles  are  taken  from  Garcia  and 
Solomon’s  (1985)  parameterization  of  mixing  induced  by  gravity- 
wave  breaking  at  45°  latitude.  The  profile  deduced  by  Allen  et 
al.  (1981)  is  also  shown  for  comparison. 
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Figure  3. 


Background  neutral  densities  adopted  in  our  calculations  for 
spring  ( - ),  summer  ( - ),  fall  (**’),  and 

winter  (-•-)  conditions  below  90  km  from  Cole  and  Kantor  (1978) 


Neutral  Chemistry 


Initialize  reaction  rates  — 
cross  sections 


Initial  guess  for  slow  species 
Noon  values  of  fast  species 


Solve  -r~  =  P.  -  L.  +  (P.  ) 
dt  1  1  1 

for  all  fast  species 

(periodic  boundary  conditions  - 

assumed  slow  species) 


Compute  diurnal  average  of  production 
loss  terms  and  densities  of  slow  species 


|  Solve  diffusion  equation  for  slow  species 


- — '  Check  con-  — » 

vergence  of  all  species  at 
_  all  altitudes - - 


lutput 


Rescale  diurnal  species  in  HO  ,  0  to  agree 

XX 

with  diurnal  average  solution  of  diffusion 
equation  -  new  guess  for  slow  species  - 
noon  values  of  fast  species 


Figure  5.  Flow  chart  for  coupled  diurnal-diffusion  program.  The  box  marked 
with  an  asterisk  (*)  corresponds  to  our  previous  diurnal  code. 

The  term  P*  represents  the  effects  of  transport  in  the  diurnal 
calculation,  as  explained  in  Rodriguez  et  al.  (1984). 


Altitude  (Km) 


MIXING  RATIO 


Figure  6.  Mixing  ratio  of  H2O  for  different  seasonal  conditions  used  in 
calculations.  The  seasonally-dependent  eddy  diffusion 
coefficients  of  Garcia  and  Solomon  (1985)  have  been  used  to 
calculate  the  H2O  mixing  ratio  for  fall  (•••)  and  winter  (-•-) 
conditions.  The  profile  calculated  by  Allen  et  al.  (1984)  (-0-), 
used  in  our  Part  I  calculations,  is  also  shown  for  comparison. 
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Figure  10.  Calculated  seasonal  variation  of  the  sodium  D  line  nightglow  at 
45°N  (histogram)  compared  to  the  observed  seasonal  variation  at 

23°S  ( - )  (Kirchoff  et  al.,  1981)  and  M°N  ( - )  (Fukuyama, 

1977).  Calculated  and  observed  emissions  have  been  normalized  to 

50R.  The  northern  hemisphere  data  of  Fukuyama  (1977)  ( - )  has 

been  shifted  by  six  months  for  comparison  with  the  southern 
seasons . 


PART  V 


Extension  of  Model  Results  to  Other  Models 


ABSTRACT 


The  diurnal  variation  and  altitude  profiles  for  Li  and  K  have  been 
calculated  using  our  photochemical  model.  Lack  of  kinetic  data  precludes 
detailed  analysis  of  the  behavior  of  the  above  species.  Our  results  are 
consistent  with  the  few  available  observations,  if  we  adopt  chemical  schemes 
similar  to  those  for  Na  and  assume  sources  equal  to  10  1  and  10  3  times  that 
of  Na. 
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1 .  INTRODUCTION 


Detailed  analysis  of  the  alkali  metals  are  difficult  because  of  the  lack 
of  observation  and  laboratory  kinetic  data.  However,  based  on  the  available 
observation  and  thermal  data,  it  is  reasonable  to  assume  that  the  source  for 
the  alkali  metals  is  also  from  meteor  ablation  and  that  the  photochemical 
sections  should  be  similar.  Thus,  a  zeroth  order  description  of  the  alkali 
metals  such  as  potassium  (K)  and  Lithium  (Li)  should  be  obtained  by  simply 
scaling  the  sodium  results  by  the  respectively  meteoric  abundances.  We  will 
discuss  some  of  the  model  findings  in  this  section  for  K  and  Li. 

2.  REVIEW  OF  OBSERVATIONS  AND  KINETIC  DATA 

The  observational  data  on  potassium  and  lithium  are  also  less  abundant 
than  for  sodium.  Twilight  observations  (Sullivan  and  Hunten,  1964;  Gault  and 
Rundle,  1968)  and  more  recent  nighttime  lidar  measurements  (Megie  et  al., 

1978;  Jegou  et  al.,  1980)  have  established  the  existence  of  K  and  Li  layers 
with  abundances  a  factor  of  10  and  103  smaller  than  that  of  Na,  respectively. 
Simultaneous  observations  of  the  Na  and  K  layers  indicate  that  the  peak  of  the 
K  layer  may  be  1-2  km  below  that  of  the  Na  layer.  The  Li  data  also  indicate  a 
peak  near  that  of  Na  (Jegou  et  al.,  1980),  although  considerable  variability 
is  observed  throughout  the  night. 

The  seasonal  variations  of  Li  have  been  observed  to  be  similar  to  that  of 
Na  by  measurements  twilight  resonance  emissions  (Gault  and  Rundle,  1968)  and 
Lidar  technique  (Jegou  et  al.,  1980).  The  observed  column  attained  a  maximum 
in  the  fall  and  winter  which  is  about  a  factor  of  five  larger  than  the  summer 
values.  It  is  difficult  to  obtain  seasonal  behaviors  of  K  because  of  the 
large  day-to-day  variations  (Gault  and  Rundle,  1968). 

Kinetic  data  for  these  species  are  essentially  nonexistent.  The  only 
available  measurements  are  those  of  Silver  et  al.  (1984)  for  the  reaction 

K  +  02  +  M  -*•  K02  +  M  ( 1  ) 

and  of  Kramer  et  al.  (1982)  for  the  Li  reaction  analogous  to  (1).  The  mea¬ 
sured  rates  (5.4  x  10  3^  (300/T)^.56  cm6  s  1  for  potassium,  3.1  x  10_30  cm^  s- ‘ 
for  lithium)  indicate  fast  three-body  recombinations  similar  to  Na. 
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On  the  basis  of  the  similarity  in  the  observed  behavior  of  the  column 
abundances,  the  chemistry  of  K  and  Li  seem  therefore  similar  to  that  of  Na. 

We  have  thus  adopted  a  chemical  scheme  similar  to  that  of  Na.  This  scheme  is 
presented  in  Table  1.  The  similarity  of  Li,  K,  and  Na  observations  suggest 
that  the  reaction 


may  occur  at  rates  similar  to  the  analogous  Na  reaction.  Thermochemical  data 
for  KO,  LiO  and  K.O2  ,  Lip2  (Murad,  1979;  Pedley  and  Marshall,  L983;  Figger  et 
al.,  1981)  indicate  that  reactions  (2)  are  exothermic.  We  thus  adopt  values 
for  (2)  ranging  from  10  to  10  cm3  s 


3.  RESULTS  AND  DLSCUSSIONS 

Calculated  profiles  for  K  at  midnight  are  illustrated  in  Figure  1.  Thes 
profiles  have  been  calculated  assuming  a  source  of  potassium  whose  magnitude 
is  one-tenth  that  of  sodium.  The  Na  source  is  scaled  to  yield  peak  densities 
of  3  x  103  cm3  s  1  for  "fail"  conditions  as  calculated  in  Part  IV.  The 
"standard"  chemistry  model  adopts  a  value  of  10  33  cm3  s  1  for  (2)  similar  to 
our  preferred  value  for  Na.  The  dashed-dot  line  and  dashed  line  denotes 
extreme  cases  with  (2)  equal  to  1<)_14  and  10_1°  cm3  s_1 .  The  Na  profile 
ral.ml.  it  ed  with  similar  chemical  scheme  is  also  illustrated  for  comparison, 
with  the  densities  scaled  down  by  a  factor  of  10.  We  also  show  a  potassium 
profile  obtained  by  averaging  our  diurnal  densities  over  the  whole  night  (thi 
solid  Line).  The  circles  correspond  to  the  data  of  Megie  et  al.  (1978). 

The  Na  and  K.  profiles  for  the  standard  cases  are  very  similar.  Agreemen 
between  our  calculated  profiles  and  the  data  is  poor  in  the  bottomside  layer. 
This  discrepancy  is  probably  due  to  the  long  integration  times  (~  12  hours) 
used  to  obtain  the  experimental  profiles.  Dynamical  processes  at  night  are 
observed  to  produce  vertical  oscillations  in  the  Na  layer  (Batista  "t  il., 
I,,3‘>).  These  oscillations  may  then  yield  average  concent  rat  ions  at  a  given 
altitude  much  higher  than  their  instantaneous  values. 


The  l, irk  ot  diurnal  data  on  potassium  precludes  deducing  any  further  con 
straiuts  on  the  unknown  reaction  rates,  part iru lar 1 v  that  of  (2).  A  similar 


situation  occurs  for  Li.  The  vertical  profiles  calculated  at  midnight  is 
shown  in  Figure  2,  again  using  the  standard  chemical  model,  a  source  equal  to 
10  ^  that  of  Na,  and  "fall”  conditions.  Variations  in  the  profile  (not  shown) 
similac  to  that  of  K  are  also  obtained  if  we  consi  I  m-  .-Kir. -me  values  for  the 
rate  of  (2).  The  shape  of  the  Li  layer  is  also  similar  to  that  of  Na  and  in 
general  agreement  with  data. 

Figures  3  and  4  show  our  calculated  diurnal  variation  for  K  (Figure  3) 
adn  Li  (Figure  4)  using  the  standard  chemistry.  These  variations,  calculated 
at  82,  86,  and  90  km,  are  similar  to  that  of  Na. 

Given  the  similarity  in  the  altitude  profiles,  it  follows  that  the  calcu¬ 
lated  seasonal  behavior  of  the  column  abundances  for  K  and  Li  are  similar  to 
that  of  the  sodium.  With  the  present  status  of  the  observational  and  kinetic 
data,  the  above  treatment  appears  to  be  an  adequate  approach.  Although  mea¬ 
surements  are  available  for  Fe+,  Mg+,  Ca+  from  rocket-borne  spectrometers,  the 
lack  of  seasonal  and  diurnal  observations  preclude  further  analysis.  Further 
clarification  awaits  more  kinetic  and  observational  data. 


Table  1:  Reaction  Scheme  and  Rate  Constants  for  Neutral  K.  Li 


M  +  O3  +  MO  +  O2 


MO  +  0  +  M  +  O2 


10  10 


lO'H 


M  +  02  +  N2  +  M02  +  N2  5.4  x  10  30 


M02  +  0  -*•  MO  +  02 


M02  +  hv  +  M  +  02 


M02  +  H  +  MOH  +  0 


MO  +  H20  MOH  +  OH 


MOH  +  hv  *  M  +  OH 


MOH  +  H  M  +  H2  0 


10  30  (T/300)  ,56<:  3.1 
10_10  -  10"13  10" 


d 

10  30 


10"10  _  10  13 


10"12 

10"10 


10"12 


Units  are  cm3  s  1 ,  except  cm6  s  1  for  (3)  and  s  1  for  (6)  and  (8). 
M  =  K,  Li. 

Measurement  by  Kramer  et  al.  (1982)  at  393°  K. 

Measured  by  Silver  et  al.  (1984). 
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Figure  Captions 


Figure  1.  Calculated  altitude  profiles  of  nighttime  potassium  and  sodium. 

The  value  of  k  (KO2  +  0)  was  set  to  10  13  ( - ),  10  14  (-••-),  and 

10  ( - )  cm3  sec  The  Na  profile  (•••)  is  calculated  using 

the  chemical  scheme  of  Table  1  and  reduced  by  a  factor  of  10  for 
comparison.  All  the  above  profiles  have  been  calculated  at 
midnight.  The  open  dot  line  (-0-)  represents  the  K  profile  with  k 

(Na02  +  0)  =  10  13  cm3  sec  1  averaged  over  the  whole  night.  We  use 

the  calculated  "fall"  background  atmosphere.  The  magnitude  of  the 
sodium  source  is  scaled  to  yield  a  peak  noontime  density  of  3  * 

103  cm  3 .  The  potassium  source  is  assumed  to  be  10  1  times  the 
sodium  source.  Circles  correspond  to  the  data  of  Megie  et  ai. 
(1978). 


Figure  2.  Calculated  altitude  profiles  of  lithium  ( - )  and  sodium  ( - )  at 

midnight.  The  value  of  k  (Li02  +  0)  is  10  13  cm3  sec  1.  The 
sodium  profile  is  calcualted  using  the  chemical  scheme  of  Table  1 
and  reduced  by  a  factor  of  1000  for  comparison.  We  use  the  "fall" 
background  atmosphere.  The  sodium  source  is  as  described  in  Figure 
1.  The  lithium  source  is  assumed  to  be  10  3  times  the  sodium 
source . 


Figure  3.  Diurnal  variation  of  potassium  species  at  (a)  82  km,  (b)  86  km,  and 
(c)  90  km.  Source  magnitude  and  background  atmosphere  are  as  in 
Figure  1.  The  value  of  K.  (KO2  +  0)  is  10  *3  cm3  sec 


Figure  4. 


Same  as  Figure  3,  but  for 


lithium  species. 
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Figure  1.  Calculated  altitude  profiles  of  nighttime  potassium  and  sodium. 

The  value  of  k  (KO2  +  0)  was  set  to  10  13  ( - ),  10  (-••-), 

and  10  10  ( - )  cm3  sec  l.  The  Na  profile  (•••)  is  calculated 

using  the  chemical  scheme  of  Table  1  and  reduced  by  a  factor  of  10 
for  comparison.  All  the  above  profiles  have  been  calculated  at 
midnight.  The  open  dot  line  (-0-)  represents  the  K  profile  with  k 
(Na02  +  0)  =  10  13  cm3  sec  1  averaged  over  the  whole  night.  We 
use  the  calculated  "fall"  background  atmosphere.  The  magnitude  of 
the  sodium  source  is  scaled  to  yield  a  peak  noontime  density  of 
3  x  103  cm”3.  The  potassium  source  is  assumed  to  be  10  1  times 
the  sodium  source.  Circles  correspond  to  the  data  of  Megie  et  al. 
(1978). 
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